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Abstract 

The emergence of a multiplet structure of the hehum-hke ions is studied within Rela- 
tivistic Schrodinger Theory (RST), a fluid-dynamic approach to the relativistic quantum 
theory of the many-particle systems. The fluid-dynamic character of RST demands to 
specify the electronic current densities for any A^-particle configuration which is ex- 
emplified here by considering the helium singlet {^Sq) and triplet (^5*1) states in great 
detail. Since the use of densities in RST is based upon the concept of wave functions, the 
new theory appears as a certain kind of (relativistic) unification of the conventional wave 
function formalism and the density functional theory, which both are the most prominent 
theoretical tools in atomic and molecular physics. As a demonstration of the practical 
usefulness of RST, the energy difference AEi\2 of the helium singlet states 2s^ and 
Is^ ^Sq is calculated for a large range of nuclear charge numbers Zex (2 < z^x < 100), 
whereas the corresponding experimental values are available only up to Zex = 42 (molyb- 
denum). The deviations of these RST results from the observational data is less than 
0,3%. 



PACS numbers: 03.65.Pm - Relativistic Wave Equations; 03. 65. Go - Solutions of Wave Equations: 
Bound States; 03.65.Sq - Semiclassical Theories and Applications; 03.75.b - Matter Waves 
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I. INTRODUCTION AND SURVEY OF RESULTS 



Surely, from the classical conception of the world, the strangest feature of the 
quantum phenomena refers to the so-called wave-particle duality . Indeed, 

such a kind of " duality" represents a certain dilemma for theoretical physics which 
always tries to erect a mathematically unique and logically coherent view upon the 
material world; whereas the notorious wave-particle duality forces us to describe 
the elementary matter sometimes in terms of fluid dynamics (wave picture) and 
sometimes in terms of point particle dynamics (particle picture) with its associ- 
ated probabilistic interpretation. The latter picture is favourable when dealing 
with statistical ensembles of particles J3|, and the wave picture is to be preferred 
for the treatment of quantum liquids j^. 

In view of such an ambiguous theoretical situation, any progress towards a uni- 
fication of both pictures should be welcome; no matter whether it refers more to 
the particle picture or rather to the wave conception of matter. Probably, greater 
efforts have been undertaken in the last decades to further develop and elaborate 
the particle picture in connection with quantum statistics, but with the recent 
experimental realization of the Bose-Einstein condensates j5| the interest in using 
fluid-dynamic concepts in quantum theory has been revived. 

Concerning the competition of both conceptual lines of thinking with respect 
to the production of testable predictions for the outcomes of real experiments, 
it should be especially interesting to consider those fields of "intersection" where 
the predictive power of both theoretical approaches may be directly opposed to 
each other and to the observational data. Naturally such a field could be atomic 
physics where experimental spectroscopy can establish very precisely the energy 
level systems of the many-electron systems. Selecting this field of physics as the 
arena of competition provides us with an additional advantage, namely the possi- 
bility to oppose the potentiality of both approaches to each other, especially when 
one has to take into account also the relativistic effects which for higher nuclear 
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charge number [z^x ^ 30, say) can no longer be neglected. Furthermore, within 
the conventional (i.e. particle) approach, one can take into account also the QED 
corrections in order to shift the predictions closer to the experimental data. This 
then must necessarily force the proponents of the fluid-dynamic picture to think 
about the incorporation of self-energy effects into their theory. Thus the arena for 
the competition seems to be well-prepared; but where are the competitors? 

On the conventional (i.e. point-particle) side, it is somewhat dissapointing to 
see that there does not exist a generally applicable and well-working relativis- 
tic quantum mechanics for many-particle systems. A possible candidate could 
have been the Bethe-Salpeter type of equations which however are afflicted with 
many deficiencies (see the critical evaluation of the Bethe-Salpeter equations in 
the literature, e.g. ref.s 0, El- Here the main problem seems to refer to the gen- 
eralization of the non-relativistic probabilistic concepts to the relativistic domain, 
which would imply to encode into the mathematical apparatus the tensor product 
of one-particle Hilbert spaces as the N-particle Hilbert space. As a consequence, 
there exist merely some auxiliary constructions in order to take account for the 
relativistic effects in atoms and molecules: Dirac-Fock approximation j^, multi- 
configuration Dirac-Fock method (MCDF) 0,^], all-order technique in relativistic 
many-body perturbation theory (MBPT) jlOl], or relativistic -^-expansion method 
For an overview of approximation methods (being based upon the tensor 
product construction), see ref. jl^ . 

However, the incorporation of the relativistic laws is no problem for a fluid- 
dynamic theory, as is demonstrated by the recently established Relativistic 
Schrodinger Theory (RST) [3, 3, which is based upon the Whitney sum 
of the one-particle bundles in place of the tensor product construction. Some pre- 
liminary tests of RST in the field of atomic physics have already demonstrated 
that the numerical predictions of this new theory are rather better than worse in 
comparison to the standard quantum theory, even if the QED self-energy effects 
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are taken into account (see ref.s 
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J3)- 

But there seems to arise a certain type of problem for the Whitney sum con- 
struction, which is not present for the conventional tensor product formalism and 
whose solution is the aim of study of the present paper: this refers to the multiplet 
structure of the many-electron configurations in atoms. In order to have a simple 
demonstration of how to proceed in more complicated situations, one may consider 
here the singlet (^5*0) and triplet {^Si) configurations of the helium-like ions with 
arbitrary nuclear charge numbers z^x (2 < z^x < 100). Neglecting first the spin- 
spin interactions and the relativistic effects, one can start with some Schrodinger 
energy level of conventional notation {ns,n's) where the principal quantum num- 
bers are denoted by n, n' and the orbital angular momentum is zero (/ = /' = 0), so 
that both (non-relativistic) electrons are in an s-state (see fig. 3 below). Switching 
now on the spin-spin interactions (and also the relativistic effects) lets the four 
different spin states of the Schrodinger configurations (ns, n's) split up into one 
singlet state (of conventional notation nis,n2S ^Sq) and three degenerate triplet 
states {nis, ^Si). This multiplet splitting of the energy levels emerges in a very 
natural way within the (non-relativistic) tensor product formalism. 

Thus, the question arises now for the new theory: In what way does this mul- 
tiplet structure emerge in RST? Observe here that the spin phenomenon in RST 
is traced back to some rotational flow of the wave function, principally not being 
different from the mechanism for the generation of the orbital angular momentum. 
Therefore the specification of the RST two-electron multiplet structure must con- 
sist in fixing the geometric pattern of the currents due to the singlet {}Sq) and 
triplet (^5*1) configurations. Our proposition for the solution of this problem is the 
following: 
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The two triplet configurations with conventional notation |S' = 1;S'2 = ±1> are 
described in RST by letting both electronic currents encircle the z-axis, either in 
the positive or negative sense. However, the "exotic" spin state jS" = 1; 5"^ = > 
is due to a rotational flow in the plane containing the z-axis so that no non-trivial 
2;-component of the total electronic spin S can arise {Sz = 0); but nevertheless 
there does exist a rotational flow {S = 1)! On the other hand, the singlet states 
^So are assumed to be due to the two currents encircling the z-axis in opposite 
directions such that both the total spin quantum number S and its 2-component 
Sz are zero: IS* = 0; 5*^ = >, see the illustration above. 

The intention of the present paper is to elaborate this proposition from the 
mathematical point of view and test its physical consequences with respect to the 
spectroscopic data. The energy eigenvalue equations for the singlet and triplet 
states are set up within the RST formalism; but the triplet case is too complicated 
to obtain a numerical solution which could be used in order to demonstrate the 
degeneracy of the three states \S = 1; Sz = 0, ±1 >. However, for the isotropic 
subset ns"^ ^Sq of the singlet states nis,n2S ^Sq one encounters a one-dimensional 
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eigenvalue problem which can be solved numerically so that a comparison to the 
experimental data is possible. More concretely, we numerically solve the RST 
eigenvalue problem for the singlet states ^Sq and Is^ ^Sq in order to calculate 
the corresponding energy difference {AEi\2) by taking the value of the RST en- 
ergy functional upon the numerical solutions. Both situations, with included and 
neglected self-interactions, are considered, and it is found that the RST predic- 
tions deviate from the corresponding experimental values ^3] by less than 0.3%, 
see table I and fig. 4 below. This result is somewhat amazing because the non- 
relativistic limit of RST coincides with the Hartree-Fock approach which misses 
the exchange effects for the singlet states (see fig. 3 below and the discussion in 
ref. Ql). Thus, the conclusion is that the exchange effect does not act in a com- 
pletely symmetric way upon the singlet and triplet states (as suggested by the 
standard theory), but is weakened for the singlet states. The mechanism of the 
RST self-interactions is discussed in detail; this elucidates the result that the in- 
clusion of the self-interactions can hardly improve the " semiclassical" value for the 
energy difference AEi/2- 

The results are elaborated by the following arrangement: 

In Sect. II a brief sketch of the fundamental RST features is presented in order 
to realize clearly the differences from the standard quantum theory. Here it should 
become obvious that RST is a rather general field theoretical framework whose in- 
dispensible features refer only to the (local) conservation of energy-momentum and 
charge, but otherwise leave great freedom to reify the general structure. Special 
realizations of RST refer to the Klein-Gordon theory of scalar particles and to the 
Dirac theory of spin-1/2 electrons; the latter realization is taken as the basis for 
the subsequent discussions. 

Sect. Ill establishes the general RST eigenvalue problem for the helium-like 
ions, see equations ()3.23j) - (j3.26dj) below. This eigenvalue problem must be com- 
plemented by the specification of an energy functional {Et, say), whose values 
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upon the solutions of the eigenvalue problem yield the desired atomic energy lev- 
els. The general structure of this RST energy functional is typical for the Whitney 
sum construction to be used for the fluid-dynamic approach in place of the ten- 
sor product construction: the total energy Et equals the sum of mass eigenvalues 
MaC^) minus those interaction contributions {^Et) which are already included 
in each one of the mass eigenvalues M^c^ and are therefore counted twice if one 
simply forms the sum of the mass eigenvalues, see equations ()3.35p - ()3.36p below. 

In Sect. IV, the general two-particle theory is cut down to the triplet config- 
urations ^Si [ortho-helium), whose RST structure is then studied in great detail. 
Naturally both states |S' = 1;S'2 = ±1> with non-zero spin-component Sz have a 
very similar geometric structure, so that it is sufficient to write explicitely down 
only one of both configurations (i.e. \S = l\Sz = +1 >). In contrast to these 
relatively simple configurations which admit a convenient non-relativistic approx- 
imation, the state IS* = 1;S'2 = > with vanishing spin component Sz has a 
distinctly different geometric structure (see the figure above) so that it seems very 
difficult to obtain an approximate solution to the corresponding eigenvalue prob- 
lem (see equations ()4.8|l - (j4.13b|l below). Therefore the question of degeneracy of 
all three states 15 = 1; = 0, ±1 > must be left unclarified for the moment (the 
degeneracy of both states IS* = 1; 5*^ = ±1 > is self-evident). 

Sect.V presents a detailed study of the singlet states ^Sq {para-helium). These 
RST counterparts of the conventional states 15* = 0; 5*2 = > have the peculiarity 
of vanishing exchange density, so that no exchange interactions (of the "electric" 
type) do exist between the two electrons. Since the "electric" interactions are 
mostly much stronger than their "magnetic" counterparts, the para-helium levels 
are not shifted (in lowest-order approximation) by the RST exchange interactions, 
which is in contrast to the predictions of the standard quantum theory being based 
upon the tensor product formalism, see fig. 3 below. Therefore it may appear as if 
RST, and the Hartree-Fock approach as its non-relativistic limit, would produce 
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inaccurate predictions of the para-helium energy levels. 

This supposition is then tested explicitely in Sect. VI by numerically solving the 
RST eigenvalue problem of Sect.V for the isotropic states ns^ ^Sq of para-helium 
for principal quantum numbers n = 1 and n = 2. This admits us to compute the 
corresponding energy difference A£'i\2 for the para-helium states 



AEi\2= E 



-E 

2^2 iSo 



(1.1) 

1^2 ISo 



and to compare our RST result to the experimental values Unfortunately 
the analogous results have been omitted by the other theoretical approaches b eing 
conveniently available in the literature, i.e. the all-order technique in MBPT [lO| 
and the -^-expansion method The result of this comparison is very favourable 
for the present fluid-dynamic approach: its predictions for the energy difference 
A£'i\2 deviate from the experimental values [Hf by less than 0.3% and do approach 
them even better for increasing nuclear charge number Zex (fig-4). This signals that 
the missing of the "electric" exchange energy for the singlet states in the RST and 
Hartree-Fock approaches may not merely be an artefact of the formalism but 
actually points to the existence of a real effect in the atomic structure! The RST 
predictions for the energy difference Ai?i\2 ()l.l|l are extended beyond the largest, 
experimentally available nuclear charge number (i.e. Zex = 42, molybdenum) where 
it must be left to future experiments to verify or falsifiate the RST results in this 
ultra-relativistic regime (Table II). 



II. RST DYNAMICS 



Similar to most of the modern gauge theories of elementary matter, the funda- 
mental equations of motion represent also in RST a coupled system of matter and 
gauge fields. In the present context, the gauge field equations are considered to be 
even more fundamental than the matter field equations since the latter type be- 
comes subjected to certain restrictions in order to guarantee the conservation laws 
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(e.g. for charge and energy-momentum); and these restrictions directly originate 
from the gauge field dynamics. Here it will turn out that the well-known wave 
equations of conventional quantum theory (i.e. of Dirac and Klein-Gordon type, 
respectively) do represent only a restricted set of possible matter equations which 
however are sufficient to ensure the desired conservation laws. Though this restric- 
tive set has an acceptable non-relativistic limit (i.e. the Hartree-Fock equations 
13, Q]), it may be nevertheless insufficient to cover all relativistic many-particle 



situations. Subsequently it is demonstrated that RST is able to provide a more 
general class of relativistic wave equations; however, for our present treatment of 
ortho- and para-helium we are satisfied with the Dirac realization of RST. 

The present procedure of generalizing the classical wave equations of Dirac and 
Klein-Gordon is based upon the use of a new field variable, i.e. the Hamiltonian 
Ti.^, which itself is a dynamical object of the theory and obeys its own field equa- 
tions {Hamiltonian dynamics) . This object Ti^ may be conceived as some kind of 
hidden variable which guides the matter through the "landscape" of gauge fields 
just in such a way that the desired conservation laws for matter are ensured. For 
the original restricted set of field equations, the Hamiltonian can be again 
eliminated, so that the conventional wave equations of Dirac and Klein-Gordon do 
reappear. However, for the more general field configurations (and especially for 
the relativistic mixtures) the Hamiltonian cannot be eliminated and remains 
an essential part of the dynamical system. 

According to the logical dominance of the gauge system, we first write down 
the corresponding dynamical equations in order to deduce thereof the restrictions 
for the motion of matter; and afterwards we introduce the Hamiltonian in such a 
way that these restrictions (i.e. conservation laws) are actually obeyed. 
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A. Gauge Field Equations 



The interactions among the material constituents of an N-fermion system are 
described by the bundle connection (i.e. "gauge potential") Afj,, which adopts its 



values in the gauge a. 
(a = l,...,iV2) 



gebra u(A^) being itself spanned by the A^^ generators Tq 



A, = -z ■ 1 + A" r„ = ('^^U^ + ^'U, . {2.1] 



Here the external potential ^'^^'^A^ refers to some external source, not being part of 
the dynamical system (e.g. fixed nucleus). For the electromagnetic interactions, 
the original gauge group U{l)xU (N) is spontaniously broken down to its maximal 
Abelian subgroup f/(l)x...xt/(l) which then describes the proper electromagnetic 
interactions, whereas the frozen gauge degrees of freedom refer to the exchange 
interactions. Correspondingly, one splits up the internal connection ^^^A^ into its 
electromagnetic part (^^"^^.4^) and exchange part (i3^), 

(^U^ = ('^"U^ + B, (2.2a) 
(-)^^ = A'^^Ta (2.2b) 
B, = B\xk - B*\xk , (2.2c) 

where the N electromagnetic generators Ta (a = 1, . . . , A^) are adopted to be anti- 
Hermitian (fa = — r^) and commuting ([ra,r;,] = 0), in contrast to the exchange 
generators Xk {k = 1, . . . , ^(^~^) Por instance, for a two-fermion system (A^ = 2) 
one has the decomposition 

A, = -I (^^U^ ■ 1 + A\n + A\t-2 + B^x - B%x. (2.3) 
with the commutation relations 

-i (n - T2) (2.4a) 

-tX (2.4b) 

IX . (2.4c) 





rs] = 


[x,x] 


[n 


X] = iX 


[n,x] 


[r2 


x] = -^x 
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Now, in order to specify the field equation for the bundle connection Af^, one 
first considers its curvature jF^j, ("field strength"), 

•^fiu ~ V^vA^ VjywA^ + (2-5) 

which, e.g., for a two-fermion system decomposes as follows 

= -^ ^'""^F^u ■ t + F'^^.T^ (2.6a) 
= -t ■ 1 + F;,r„ + G,,x - G\^x , (2.6b) 

with the curvature components being given in terms of the connection components 
as follows 

^'""^u = (^^'U, - V. (^^'U^ (2.7a) 

= V,A\ - V,A\ + z {B,B: - B,B*^) (2.7b) 

= V,A\ - V,A\ - I {B,B: - B,B%) (2.7c) 

G,, = V,B, - V^B, + I {A\ - A\)B, - i {A\ - A\)B, (2.7d) 

= V,B: - V,B\ - I {A\ - A\)B: + t {A\ - A\)B\ . (2.7e) 

Once the curvature jF^j, of the connection A^ has been introduced, the field equa- 
tion for Afj, is selected to be the usual (non-Abelian) Maxwell equation: 

'D'^T^,, = -AmasJ, (2.8) 

{VxT^u = VxJ^^u + [Ax,J^f,u]) ■ 

Decomposing here the (Hermitian) current operator J'^j, with respect to the chosen 
Lie algebra basis {tq,} as 

J, = ■ 1 + zr^T^ , (2.9) 
i.e. more concretely for our two-fermion system 

J, = ^^"^Jm ■ 1 + ^JV^ + Kh*^X - Kx) , (2.10) 
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one finds the component version of tlie operator-valued RST-Maxwell equations 
()2.8|) to be of the following form: 









(2.11a) 




- B'^^G^,) 


= ^T^asj^y 


(2.11b) 




- B*^'G^,) 


= Anasfy 


(2.11c) 


V^G,^ + z {A\ - A\)G^, - I (FV 


- K>^)B^ 


= An as fy 


(2. lid) 


V''G%, - I {A\ - A\)G*\ + I {F\, - 


- F\,)B*' 


= -An as fy . 


(2.11e) 



At this stage of the development of the theory, one necessarily becomes con- 
fronted with the first restriction upon the motion of matter which itself is thought 
to produce in some way those Maxwell currents j"^ ()2.9|) which enter the Maxwell 
equations ()2.8p . The point here is namely that any bundle curvature JF^^ ()2.5|) 
must obey the identity 

[V^,V,]J^'''' = (2.12) 

which by adopting that link ()2.8|) of the Maxwell equations to the current J'f^ 
immediately implies the "continuity equation" in operator form 

= . (2.13) 

Actually this is a severe restriction upon the motion of matter which becomes 
more obvious when written in component form, e.g. for the present two-fermion 
system: 

V''j\ = z{B^j\ + B*^f^) (2.14a) 
V^f^ = -^(B'^J\ + B*^J\) (2.14b) 
3\ = -I {A\ - A\)j'^ + I B^ {j\ - j\) (2.14c) 
j\ = ^ {A\ - A\)f^ + z B*^ {j\ - (2.14d) 

This source system is to be interpreted in the following way: for given Maxwell 
currents j"^ with sources (V'^j'^^) as specified by the present equations ()2.14aj) - 
()2.14d|) . the connection components {A°-^; S^} must be thought to be solutions of 
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Fig. 1: Self- Consistency of Currents and Potentials 



the RST-Maxwell system (j2.11bj) -( )2.11e|l with just these currents j"^ emerging on 
the right-hand sides! (see fig.l). 

A similar self-consistency problem for the coupled matter and gauge field dy- 
namics arises for the energy-momentum density of matter {^^^T^^, say): According 
to the general logical structure of the gauge field theories, one expects that the 
well-known Lorentz force density ^^^fu should act as the generator of the matter 
energy-momentum, i.e. 

V'^(^)T^. = (% . (2.15) 

For the present situation of a composite system (N particles), the Lorentz force 
will consist of an external part ^^^^f,^ and an internal part ^^^fu 



(2.16) 



Here the external part ^^^•'/^ is assumed to be built up by the external field 
strength ^'^^^F^^ and the total electromagnetic current as the sum of the electro- 
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Fig. 2: Energy- Momentum Balance 
magnetic gauge currents = 1, A^) 

N 

i.e. one puts 



a=l 



A similar construction suggests itself for the internal part '•'^^/^ 



(2.17) 



(2.18) 



(2.19) 



which however obliges us to introduce a covariantly constant fibre metric for 
the Lie algebra bundle in order to lower and lift the gauge algebra indices 



3 an 



Kap ff, etc. 



(2.20a) 
(2.20b) 



But here it should be self-evident, that the requirement of energy- momentum bal- 
ance ()2.15|) poses again the same logical problem as it already did emerge from the 
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source equations for the currents ()2.14a|) - ()2.14d|) : namely to let the matter move in 
such a specific way that the source '■'^-*/^ of its energy-momentum content '•^^-T^jy be- 
comes the product of its currents Jq,^ and the field strengths being generated 
by just those currents ja^ via the RST-Maxwell equations (j2.11bj) - ()2.11ej) !(fig.2) 

A very pleasant feature of RST is now that both self-consistency problems ()2.14|1 
and ()2.15|) for the motion of matter can be solved in one step, namely by setting 
up the right equation of motion together with a suitable definition of currents j"^ 
and energy- momentum The interesting point here is that the conventional 

field equations of matter (i.e. Dirac and Klein-Gordon) do reappear only as a 
restricted subset of a larger set of possibilities. 

B. Matter Dynamics 

The RST description of matter is based upon the intensity matrix X, which is 
a Hermitian (4A^ x 4A^)-matrix (X = X) for a system of fermions. If this system 
is to be described by a pure state the intensity matrix X degenerates to the 
tensor product of 

X^*®^. (2.21) 
The field equation for X is the Relativistic von Neumann Equation (RNE) 

v^i=^jin,,-n,,i] (2.22) 

or for the situation of a pure state \Ef, respectively, the Relativistic Schrodinger 
Equation (RSE) 

ihcV^-^ = n^-^ (2.23) 

The intensity matrix X(x) and the wave function \E^(x) are conceived as sections 
of appropriate fiber bundles with the "potential" ()2.1|) working as the bundle 
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connection. Therefore the first task is to ensure the generally valid bundle identities 



[V^V, - V,V^] 1 = [J^^,, 1] (2.24a) 
[V^V, - V,V^]^ = T^,^ (2.24b) 

by an appropriate choice of the field equations for the Hamiltonian (Eamilto- 
nian dynamics) . Consequently, the first field equation for the Hamiltonian is the 
integrability condition 

which in combination with the matter dynamics ()2.2H) - ()2.22|1 just ensures the 
validity of those bundle identities (j2.24a|) - (l2.24bj) . 



C. Charge Conservation 

The second field equation for the Hamiltonian 7i^, which is a g[(4A^, C)-valued 
one-form, must now refer to the guarantee of the conservation laws discussed in 
the preceding subsection. For this purpose, one first introduces the total velocity 
operator as a Hermitian one-form (F^ = F^) which takes its values also in the 
Lie algebra gl(4A^, C). By means of this object, one alternatively defines the total 
electromagnetic current ()2.17|) through 

j^ = tr{l-T^) (2.26) 

and then one requires this current to obey the conservation law 

V^j^ = 0. (2.27) 

Indeed, adding up both source equations ()2.14ap and ()2.14b|) for the two-particle 
situation (A^ = 2) actually yields a result which just coincides with the present 
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conservation law ()2.27|) . Thus, in this special respect, our choice of the matter 
dynamics would be consistent with the gauge field dynamics, provided we suceed 
in deducing the conservation law ()2.27|1 also from the proposed matter dynamics 
()2.22|1 . However, this can easily be done by simply transcribing the derivative 
(|T77jl of the current (jT^ to the intensity matrix I and using the RNE fTI^ 
(or the RSE ()2.23|) for the case of a pure state) which ultimately yields 

= V^j^ = tr I J ■ (v'^T^ + [n^T^ - T^n^] ^ I . (2.28) 

Thus, in order to ensure the validity of the desired conservation law ()2.27|1 for 
any matter arrangement X, we require the total velocity operator to obey the 
following operator equation 

'D'^T^ + [H/^r^ - T^n''] = . (2.29) 

This equation is to be understood in the sense that for given Hamiltonian 7i^(x) 
(as solution of the Hamiltonian dynamics) the total velocity operator must be 
a solution of its source equation ()2.29j) . 

Next, one turns to the operator- valued continuity equation ()2.13|) which reads 
in component form 

D^ja>. = , (2.30) 

with its contravariant two-particle realization being given by equations ()2.14ajl - 
(|2.14djl . Here the covariant derivative (D) of the gauge objects is defined through 



Df^jau = '^f.jau - ^^^a^iinu (2.31a) 
D,j\ = V,r, + uj"p,f,. (2.31b) 

where the connection one-form = |ci;"g^} takes its values in the adjoint repre- 
sentation of the gauge algebra, i.e. 

o;^ = C^,aA\ , (2.32) 
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and the structure constants C*"/3-y are defined through the usual commutation re- 
lations 

W.M=C\pT,. (2.33) 

Now in order to ensure the desired validity of the source equation ()2.30j) for the 
gauge currents Jq,^, one introduces the gauge velocity operators Va^ which then 
build up the currents ja^ through 

jaf, = tr{X ■ Va^j,) , (2.34) 

similarly to the construction of the total current ()2.26p . Accordingly, the tran- 
scription of the derivative (D) of ja^ in equation ()2.30p to the gauge velocity 
operators Va^i ()2.34j) yields a result quite analogous to equation ()2.28|1 . i.e. 

Therefore, in order to satisfy the requirement (j2.30|) . one adopts the velocity op- 
erators to be solutions of the following source equation: 

V^v^^ + ^af. - v^^ ■ W) = 0, (2.36) 

where the covariant derivative of the gauge operators is defined through 

T>\Va^, = VxVaf, + [Ax, Vaf,] " ^'^a\^l3^i ■ (2.37) 

However, in order to not introduce too many new objects requiring their own 
field equations, one would like to express those gauge velocity operators v^^ in 
terms of the already existing objects Tq, and F^, i.e. we try to identify v^^ as the 
anticommutator of the gauge algebra generators Tq- and the total velocity F^: 

Va^.= '-{raJ>.] . (2.38) 

Since the electromagnetic generators (a = 1,...,A^) are adopted to be anti- 
Hermitian (f^ = — t^) and the velocity operator is taken to be Hermitian (F^ = 
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r^), the electromagnetic velocity operators (a = 1, . . . , A^) are also Hermitian 
{va^ = Vap) and consequently the electromagnetic currents ja^ {a = 1,...,N) 
()2.34|1 are real, just as the electromagnetic potentials A"'^ and their field strengths 
F"-^j^. In contrast to this, the remaining N(N-l) exchange generators {Xk,Xk', k = 
1, . . . , ) are arranged in ^^^^"^^ pairs of Hermitian conjugate operators, so 

that the exchange potentials B^^ and their field strengths G^^,^ are complex, as 
well as the corresponding exchange currents h^^; see for instance the two-particle 
realization ()2.11b|) - ()2.11ep of such an arrangement. 

Introducing the chosen ansatz ()2.38j) for the velocity operators Vafi into their 
field equation ()2.Hfi|l must then necessarily yield some condition upon the total 
velocity operator F^, namely one finds the following combination of commutators 
[ .. , .. ] and anti-commutators { .. , .. } : 

^ {r„, OT, + ^ {mr, - T,n^) I + ^ [r„, mT, + T,W] - (2.39) 

Thus, summarizing the situation with the charge conservation laws, it becomes 
clear that both conservation laws ()2.27j) and ()2.3()|1 can simultaneously be satisfied 
by requiring the total velocity operator to obey the former equation ()2.29|) and 
let it commute with the gauge algebra generators Tq, 

[r«,F^] = 0. (2.40) 

However, observe here that these two conditions are actually sufficient to guarantee 
the total charge conservation (j2.27|) but leave us with one remaining term for the 
gauge continuity equation ()2.30|1 . More concretely, introducing the (Hermitian) 
mass operator Ai (= Ai) through 

Hf'Tp + T^n^ = 2 Mc^ , (2.41) 
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then the source equation for the gauge velocity operators v^^ becomes in place of 
the desired result ()2.36p : 



n^] [ra,M] 



(2.42) 



and consequently the source equation ()2.30j) for the gauge currents ja^ acquires an 
inhomogeneous term: 

Mc' 



D^jat. = tr\l 



h 



(2.43) 



Clearly, one can here easily bring down to zero the right-hand side in order to attain 
the desired homogeneous equation, namely by simply letting the mass operator M. 
being proportional to unity 

. (2.44) 



D. Energy-Momentum Conservation 

Before one can make a final decision for the choice of the total velocity operator 
as a solution of its source equation ()2.29j) . one must first specify the second 
field equation {conservation equation) for the Hamiltonian Ti^. Here one expects 
that this equation for Ti^ will emerge in connection with an energy-momentum 
conservation law for matter. More concretely, one tries to introduce the energy- 
momentum density of matter {^^'^^^y, say) in such a way that its source is just 
the well-known Lorentz force density '^^^fy ()2.1(j|) due to the external and internal 
gauge fields. This idea may be realized by first defining an energy-momentum 
operator for fermionic matter (^*^^j/, say) through 

(^^>r^, = i {t;h, + n,T^ + t,h, + n^v,] , (2.45) 

and then the energy-momentum density for matter ^^'^^T^u is obtained quite simi- 
larly as for the current densities ()2.2(i|l or ()2.34j) by putting 

(*^>r,. = tr (J ■ (^■^>r,.) . (2.46) 
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Whether (or not) this is a successful construction is decided by looking for the 
divergence of this tensor: 

^'%. = tr I J (v^ ^'%^ + ^ [W ^'%, - ^1 ) } . (2.47) 

As mentioned above, cf. p.l5|l . this divergence of the matter tensor ^^'^^^u should 
equal the Lorentz force density ^^^fu which surely will be true when the " RST diver- 
gence" of the operator ^^^^^u yields just the corresponding Lorentz force operator 
(^)f., i.e. 

^^%. + [W ^'%. - = (^)f. , (2.48) 

so that the desired force density fu ()2.1(i|l emerges in the usual way as 

(% = tr(j-(^)f.) . (2.49) 

However, the desired RST divergence of the operator can be computed 

in a straight-forward manner under use of the integrability condition ()2.25j) and is 
then ultimately found to look as follows: 

^'%. + [n^ - ^'%.n^] (2.50) 
= ^ jz^c (r^ j-^, + ^^,r^) + (v^n,, - n^n^^ + (vn^ + n'^n^^ 

+ [v^v, - v,v^] -w + w- [v^T, - v,T^] + V, [v^n^ - TY^r^] | . 

Here, according to the original idea ()2.48p . the right-hand side must yield the 
Lorentz- force operator ^^^f^, which is to be conceived as the product of velocity 
operators and curvature components, see the first term. Obviously, in order to 
achieve this goal, one can resort to the previous equation (j2.29j) which puts the 
RST divergence of the total velocity operator to zero; but this is of course not 
sufficient, and further assumptions concerning and must be adopted. Any 
such choice, which then validates the charge conservation laws ()2.27|) . ()2.30|) and 
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the energy- momentum conservation (or balance, resp.) law ()2.15|) . is considered as 
a particular realization of RST. We now turn to the most important one of these 
realizations. 

E. Dirac Realization of RST 

Evidently the requirement of conservation laws can be expressed in terms of 
certain conditions upon the velocity operator and the Hamiltonian First 
consider the case of F^ which must have vanishing RST divergence ()2.29j) and is 
required to commute with the gauge algebra generators r^, cf. ()2.40|) . in order that 
the charge conservation laws be satisfied. An additional condition arises from the 
requirement of energy-momentum balance ()2.48|) and (j2.50p where we assume now 
also the vanishing of the curl of F^ 

Vf,T, - V,T^ = . (2.51) 

However, we even want to go one step further and put the covariant derivative of 
F^ to zero 

V^T, = . (2.52) 

Fortunately this implies a great simplification of the mathematical formalism, be- 
cause, on account of the commutation relations ()2.40j) . one can always choose a 
gauge where (at least locally) the velocity operator F^ is absolutely constant: 

d^T, = . (2.53) 

Thus, the velocity operator becomes a fixed object, independent of any field kine- 
matics and dynamics. For the Dirac realization of RST, one chooses for F^ the 
A'^-fold direct sum of the well-known Dirac matrices 7^ i.e. 

T;, = 7m ® 7m ® ■ ■ ■ ® 7m • (2.54) 
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Naturally, such a choice of the velocity operator must imply further conse- 
quences. Returning for the moment to the vanishing RST divergence ()2.29p of F^, 
one deduces from that equation the following relationship between and H^: 

n^'-T^-T''-n^ = . (2.55) 

Effectively this is a condition of Hermiticity and thus one can simplify the former 
mass relation ()2.4H1 to 

H'^T^ = V^Hf, = Mc^ . (2.56) 

Furthermore, for the Dirac realization of RST for a system of identical particles, cf. 
()2.44|) . the proper gauge continuity equation ()2.30|) is implied by the more general 
relation ()2.43|) . For non-identical particles, one puts the currents and potentials 
of the exchange type to zero, so that the source system ()2.30|) becomes reduced to 
N true conservation laws (a = 1, . . . , A^) 

V^f^ = . (2.57) 

Indeed, this is obtained from the general relation ()2.43j) by linear combination of 
the mass operator M. from the electromagnetic generators alone: 

M = i M^Ta , (2.58) 

with real- valued particle masses Ma- In the present paper, we however treat exclu- 
sively identical particles and will therefore always apply the identical-mass relation 
()2.44|) in combination with the Dirac realization. 

Finally, the source equation of the Hamiltonian Ti^ must be determined in such 
a way that the right-hand side of the divergence relation ()2.50|1 for the energy- 
momentum operator ^^'^^^u can really adopt the desired Lorentz form as indicated 
by the former claim ()2.48j) . However, this goal can be attained by adopting the 
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following source equation for ( conservation equation ) 

v^n, - ■ n, = -^/icj (^^^ + s^'^j-^.j (2.59) 



Indeed, by this choice of the source of Ti^ the RST divergence of the energy- 
momentum operator ()2.50|) becomes 

^'%u + ;^ ■ ^'X'^ - ■ n^) (2.60) 
= ^ (^^" [s^A, r,] + {T,,, rn) = z^c^^.r'^ , 

where the following commutation relation of total velocity and the Spin(l,3) 
generators T.^^, is to be observed: 

[S^A, r^,] = Tf,gxu - Txg^,u ■ (2.61) 

Thus with this result the Lorentz force operator ^^^f^ ()2.48p is found to be just of 
the expected form, namely as the product of field strength jF^j, and velocity 



^, = ihcT^,-T^ . (2.62) 



This pleasant result for the force operator ^^^^u transcribes then immediately to 
the corresponding force density ^^^fv ()2.49j) 

= (^"V. + ^'^U , (2.63) 

namely by simply inserting the curvature decomposition ()2.6a|) and observing the 
trace definitions for the total current fl2.26|) and the gauge currents j^^ fl2.34|) . 
Thus the former assertions ()2.18j) - ()2.19j) are actually validated. 

For the special situation where matter is in a pure state, cf. fl2.21|l . the RSE 
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()2.23|) may be converted to the well-known Dirac equation for the N-particle wave 
function \1/ 

ih TO^* = Mc-^ . (2.64) 

For this purpose, one merely has to differentiate once more the RSE ()2.23j) under 
use of the original form of the conservation equation ()2.56|1 . The alternative version 
fl2.59|l of the conservation equation can easily be deduced from its original form 
fl2.56|) and is needed when one wishes to transcribe the Dirac equation ()2.64|) to 
its Klein-Gordon counterpart as a second-order wave equation: 

-Y) ^ = - S^'"-^^,^ . (2.65) 

For the subsequent applications of the theory to the two-particle systems, it is 
useful to write down the Dirac equation ()2.64j) also in component form: 

ih YD^ipi = Mc^i (2.66a) 
thYD^iJ2 = Mcij2 , (2.66b) 

where the gauge-covariant derivatives of the single-particle wave functions ipa (a. = 
1,2) are given by 

^mV'i = d.^P, - I [(^^U^ + A\] ^1 - I B^^2 (2.67a) 
D^^Ij2 = d^i}2 - t [('^"U^ + A\] - I B\^^ . (2.67b) 

This coupled Dirac system will readily be specialized down to its stationary form 
in order to set up the energy eigenvalue problem for the two-electron atoms. 



III. ENERGY EIGENVALUE PROBLEM 



In order to convince oneself that RST can predict qualitatively the correct mul- 
tiplet structure of the atomic spectra, one will first consider a very simple situation, 
i.e. the helium problem (two electrons around the fixed nucleus). Here one may 
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further simplify matters by treating only those stationary two-particle states which 
have vanishing angular momentum (L = 0), i.e. the singlet and triplet states ^Sq 
and ■^Si. 

The standard (non-relativistic) quantum mechanics (in lowest order of approx- 
imation) says that the four possible combinations of the two spins 

|T)®|T), U)®U), |T)®|i), U)®|T) 

will be arranged in such a way that there arises a singlet state 

l^/) = ^{|T)®li)-|i)®IT)} (3.1) 

which has vanishing total spin quantum numbers {S,Sz) = (0,0), being defined 
through 

S^l^Ji) = SiS + l)fi'\i;i) (3.2a) 

S,\^i) = S,h\^j) (3.2b) 
{-S <S,<S) , 

and furthermore there arise three triplet states ) with = 1 and = ±1,0 



5, = +l: |<)) = |t)®|T) (3.3a) 

^. = 0: I<^) = ±{|T)®U) + U)®IT)} (3.3b) 

S. = -l: lV'fe'^) = li)®U). (3.3c) 

Since for these two different kinds of states and \ipni) the single-particle spins 
S{a) (a = 1, 2) do build up the total spin S in two different ways (namely "parallel" 
for the triplet states S = 1, and "anti-parallel" for the singlet state S" = 0), it 
does not come as a surprise that the electronic interaction energies of the two spin 
arrangements are different: the energy eigenvalue of the singlet states (para- 
helium) is increased (in lowest perturbation order) by the "exchange energy" AEc- 

27 



Ai?(-') = 11 d\-£f' ^-^'^^^"[^^y^'^^^"^'^^^ , (3.4) 

whereas the energy of the triplet states \ipin) (ortho-hehum) is loweredhj the same 
amount a4"'"'\ Here the raising and lowering of energy refers to the "unper- 
turbed" levels where the spin effect is completely neglected and the corresponding 
unperturbed wave functions '?/^„(f^ refer to the (spinless) Schrodinger eigenvalue 
problem. Consequently the ortho/para level splitting according to the standard 
theory 6stE becomes twice the exchange energy 

SstE = Em -Ej = -2 a4"'"') , (3.5) 

which is in best agreement with the experimental situation (see any standard 
textbook [l| and also the discussion in ref. |3]). 

Concerning now the situation in RST, which has the well-known Hartree-Fock 
approach as its non-relativistic limit, it is well-known that the latter approximative 
method does not always produce useful predictions, though being of great help in 
many respects (see the discussion of this point in ref. j20|). For the present case 
of ortho/para splitting 6E the ordinary Hartree-Fock method (and therefore also 
RST) predicts only half of the correct value SstE ()3.5|1 

SHFE=^6stE = -AE^^''''^ . (3.6) 

The reason is here that the ordinary HF method yields an exchange energy AE'c"'" 
fl3.4|l exclusively for the triplet states \ipin) (ortho-helium), but not for the singlet 
states see fig.3. Thus the HF method appears to be helpful for computing 
the ortho- level system but not for the para- level system of helium. Clearly, since 
RST is the relativistic generalization of HF, this methodological ambiguity must 
transcribe to RST which thus could be expected to be a useful technique only 
when dealing with the ortho-states ^Si. However the situation becomes more 
favourable also for the para-states "^Sq when one considers high-energy situations 
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(i.e. high nuclear charges Z » 1) because in this highly relativistic regime it is 
more important to take into account the relativistic effects than the correlation 
effects beyond Hartree-Fock (see below and ref. for a demonstration of this). 

The ortho/para energy shift comes about mainly through separation/crowding 
of the electronic charge clouds by (anti-) symmetrizing the spatial parts of the 
wave functions relative to the symmetry properties of their spin parts, according 
to the Pauli antisymmetrization postulate for the total wave function. 

In any case it is necessary and desirable to give an explicit demonstration of 
how the singlet and triplet structure of the helium-like atoms does arise within 
the framework of RST, so that the problem of the missing exchange energy for the 
singlet states can be elucidated, see the discussion below equation ()4.18j) . 
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Standard theory 



HF, RST 



6hfE = ^6stE 



-1 = 



ns n s 



=+1 



SstE 



Fig. 3 Ortho/Para Level Splitting 

The standard level splitting 6stE (|3.5j) of a Schrodinger state ns,n's is twice the 
Hartree-Fock value ShfE because the latter approach treats electrons with 
different spins 1 1 ) | J, ) (para-helium) as non-identical particles missing the exchange 

effect. 



A. Stationary Two-Particle States 

For a stationary state, one expects that certain physical observables (e.g. en- 
ergy) are time-independent, whereas other objects of the theory are admitted to 
undergo regular oscillations. Thus the two components ipa{^,t) (a = 1,2) of the 
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wave function \1/ are written as 

V'i(r,t) = exp 

'ilj2{r,t) = exp 



.M2C^ " 



h 



(3.7a) 
(3.7b) 



where the time-independent Dirac spinors ipai"^) may be further decomposed into 
two-component Pauh spinors '^'^^(j)±{r) according to (a = 1,2) 



^a(^) 



(3.8) 



Furthermore, each of the four Pauh spinors '^'^)0-|-(r) may be decomposed with re- 



spect to an orthonormal basis, namely either the basis C,^ for the upper components 



or, resp., the basis C,^^ for the lower components 



(«), 



.(r) = -I ^"^R. ■ Cl+^ - I • Cj"^ 



(3.9) 



(3.10) 



iere, the two basis systems are defined in terms of the usual standard basis ( ^ 



3,rn 



2l| as 



1 1 

'2 





(3.11a) 
(3.11b) 



1,1 1,1 



Recall here, that the single-particle standard basis ^ , C^l ^ has the well-known 
eigenvalue properties 



J2 ^3,m 


= jU + i)h\T 


(3.12a) 


J. CT 


= mh C'^'i"^ 


(3.12b) 




= i{i + i)fe c'r 


(3.12c) 




= lil + i)^'CT- 


(3.12d) 
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According to such a stationary ansatz for the wave functions, one expects that 
also the corresponding gauge fields are either time-independent or regularly oscil- 
lating: 



B^{f,t) =^ exp 



.M1-M2 2, 
2 : C t 



■ \B.{r).-B{r', 



(3.13a) 
(3.13b) 



Moreover, the RST-Maxwell equations (j2.11bj) - ()2.11e|l clearly demonstrate that the 
currents are intimately linked to the gauge potentials and therefore they may be 
expected to be both of a similar form. Thus one assumes that the electromagnetic 
currents j""^ are also time independent: 



fa = -Kfi = i^a 7m i^a 



(%(r-) 



, -kair) 



(3.14) 



whereas the exchange current /i^ ( = js^j = — J^^) oscillates countercurrently with 
respect to the exchange potential fl3.13b|) : 

' M1-M2 2 / 



hf, = ipi 7m ^2 



exp 



h 



■cH 



ho{r) , -h{r) 



(3.15) 



Concerning the external potential 1)2.11) . we shall treat the helium problem 

where ^^^^A^ is due to the Coulomb potential of the nucleus 



{(^^Uo(f),0} 



(3.16a) 
(3.16b) 



he 



According to the restriction to such a spherically symmetric potential '^^^'•'A^, one 
will resort to the use of spherical polar coordinates (r, ip) and will decompose 
all three-vectors with respect to such a moving basis system {er,e^,e^}, i.e. one 
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puts 




(3.17d) 



(3.17b) 



(3.17a) 



(3.17c) 



B. Currents 

The current densities j"^ = {ka^; h^} are especially interesting as far as the 
geometric difference between an ortho-field and a para-field configuration is con- 
cerned. Indeed one has to expect that the physical difference between both field 
configurations is ultimately a consequence of the different geometric patterns of 
the associated ortho- and para-streamlines. The point here is that the spin in RST 
is traced back to the rotational flow generated by the wave functions and therefore 
the spin-spin interaction emerges as the magnetic field energy of those different 
rotational flows. 

Inserting the general ansatz ()3.9|) - ()3.10|) into the RST currents ja^ ()2.34|) 



■2 



k2ti = ^27m^2 



(3.18a) 




kit, = V^i7mV^i 



(3.18b) 



J3m 




(3.18c) 



J4M 




(3.18d) 



yields for the electrostatic charge densities *^°^fco(^ ()3.14p (a = 1,2) 
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and similarly for the corresponding electric three-currents ka{r) ()3.17ap (a = 1,2) 

(<')k^ = ^{^"^Rl ■ ^"^R- - ■ ^"^R*^ + ^"^Sl ■ ^"^S^ - ^"^8+ ■ ^^^S*J] (3.20a) 
= ^{^"^Sl ■ ^"^S. + ^"^8+ ■ ^"^S*, - ■ - ■ (-^^i?!} (3.20b) 

(3.20c) 

Furthermore the exchange density ho{f) becomes 

h, = ^{^'^Rl ■ ^'^R+ + ^'^Sl ■ (2)5+ + (^)i?*„ ■ + (1)5* ■ (2)5_} (3.21) 

and the exchange current h{r) (|3.17b|l is found as 

= ^{(1)/?; ■ + (1)5; ■ (2)5_ - . (2)/?^ _ (1)5* . (2)^^} (3.22a) 
= J-{WSI . + (1)5: ■ (2)5+ - . _ Wr*_ . (^)R^} (3.22b) 

= ±{e'YK ■ '^'^S- + (i)/21 ■ (2)5+] - e-'T'^Sl ■ ^'^R. + (^^^^ ■ ^''>R+]} . 

(3.22c) 

These densities and currents are now relevant for the energy eigenvalue problem 
insofar as they determine the potentials A"^ via the solutions of the RST Maxwell 
equations (j2.11b|l - ()2.11e|l : these potentials then enter the covariant derivatives 
D^ipa of the wave functions ()2.67ap - ()2.67b|) which in turn enter the eigenvalue 
equations ()2.66aj) - ()2.66b|) . This will readily become more evident when writing 
down now the stationary form of the coupled Dirac system ()2.66ap - ()2.66b|) . 

C. Mass Eigenvalue Equations 

Naturally, the mass eigenvalues Ma fl3.7a|) - ()3.7b|) will form an essential con- 
stituent of the atomic energy levels {Et-, say), but they are not directly identical 
with those atomic energies. Rather, one will have to specify the atomic energies 
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Et in terms of an appropriate energy functional, which, besides by the mass eigen- 
values (MaC^), is built up also by the energy content of the gauge fields. However, 
the energy contribution of the matter fields is essentially determined by the sum 
of mass-energies (Mic^ -|- M2C^); and therefore the mass eigenvalue problem must 
first be solved before the value of the energy functional Et upon the corresponding 
solutions of the coupled matter and gauge-field system can be calculated. 

Once the stationary form of the wave functions has been specified through 
equations ()3.7|) - ()3.12|) . one simply inserts this stationary ansatz into the coupled 
Dirac system ()2.66a|l - (j2.66bj) and then finds the corresponding mass eigenvalue 
equations in terms of the wave amplitudes ^°'^R± and ^"'^S±. These equations may 
perhaps look a little bit lengthy when written down in spherical polar coordinates, 
but it is worthwile to display them explicitly because a lot of information about 
the peculiarities of the electromagnetic and exchange interactions can be drawn 
from them. The equation for the wave amplitude reads 



am. 



% 



dr r 



e 



d-t} 2 sin 



e 



r sin-i? 



i - 1 ■ ^'^s. 



dip 2 



(3.23) 



- te''> ■ ■ (1)5+ + [B^ + zi?,] ■ (2)/?^ - ^e^^B^ ■ ^'^S^ = _Ml±^c ■ ^'^R_ . 

n 

Obviously, the "positive-energy" component ^^^R+ of the first particle's wave func- 
tion couples via the electrostatic potential ^^^Aq and electric exchange potential 
Bq to the "negative-energy" components ^^^R^ and (^■'-R_, respectively; whereas 
the coupling to the positive-energy components ^^^R+ and (^)i?+ occurs via the 
magnetic potentials ^"^^A^ and (^^A,,. In contrast to this, the coupling of to 
the " secondary" positive-energy components ^'^^S+ occurs always via the azimuthal 
components ^"^^A^ and B^, respectively. Of course, the coupling of the negative- 
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energy component follows an analogous pattern: 



dr 



d e 
+ 



2 sin-i? 



+ 



r sin-i? 



2 



(3.24) 



Thus it becomes evident that the interactions of the electric type couple the wave 
amplitudes of different energy types, whereas the magnetic interactions connect 
the amplitudes of the same energy type. 

Another striking feature of these two-particle interactions concerns the coupling 
of the "primary" amplitudes ^"'^R± to their "secondary" counterparts ^'^^S±, and 
vice versa. In order to see this more clearly, one writes down also the eigenvalue 
equations for the "secondary" components ^^^S±: 



dr 



dd 



-i-d 



2 sin?? 



-i-a 



r sin-i? 



9 ^'K ^ (Dp 



+ [(^^Uo + ('Uo] ■ (i)S_ + Bo ■ - - 2 ■ (2)^,] . (1)5, 



+ le-'^ ■ (2)A^ • (i)i?+ - [B^ - iBr] ■ (2)5, 



le 



-id 



B.„ . (2)/2, 



Mi + M 



(3.25a) 



• (i)5_ 



d (i)5_ 
dr 



- ■ (i)5„ 
r 



d (1)5, 



2 sin?? 



■ (i)5_ 



r sin-i? 



d (1)/?, 

dip 



- ■ (i)i?_ 
2 



[(^^■)Ao + (')Ao] ■ (1)5+ - 5o ■ + [(2)A^ + I ■ (2)^,] ■ (i)5_ 



(3.25b) 



le 



-i-d 



■ (2)A^ ■ (i)/2„ + [5^ + zBr] ■ (2)5_ 



Ml -M 



■ (1)5. 



From here it is obvious again that the electric potentials {(")Ao, Bq} mediate 
between the positive/negative energy components ((")5+ ^ (")5_) whereas the 
magnetic potentials {Aa, B} mediate between the components of the same energy 
type ((")5+ ^ (")5+, etc.). However, what may also be seen clearly from the 
complete set of eigenvalue equations for the first particle (a = 1) ()3.23|) - ()3.25b|) 
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is the fact that the couphng of the primary (R) and secondary (S) components 
{R ^ S) runs along the azimuthal direction. This will readily become important 
when we subdivide the whole set of solutions for the eigenvalue equations into the 
subsets of ortho-type and para-type, respectively. Finally, let us mention that, for 
the sake of brevity, we will not write down here the eigenvalue equations for the 
second particle (a = 2), because this is obtained from the present first-particle's 
equations ()3.23p - ()3.25b|) simply by means of the following replacements: 

Ml ^ Ma (3.26a) 
(1)/?^ ^ (2)/?^^ (1)5^ ^ (2)5^ (3_26b) 

(^Uo ^ ^^Uo, Bq ^ B* (3.26c) 
Ai ^ A2, B ^ B* . (3.26d) 

Two remarks must be made concerning the solutions of the mass eigenvalue 
problem ()3.23|) - ()3.26d|) . Namely first, these solutions are unique only if one imposes 
some normalization condition upon the single-particle wave functions ipai'f^, i-e. 
one needs some normalization condition upon the wave amplitudes ^"'^R±{r) and 
^"■'>S±{f). The desired condition may be deduced in a very natural way from the non- 
abelian RST-Maxwell equations ()2.11b|) - ()2.11ep : Recasting the electromagnetic 
part ()2.11b|) - ()2.11c|) of this system into the following form [l5| 

d^F\, = -47ra, ^% = -Arras + ^ [B^G%, - fi^^G^.] | (3.27a) 

= -47ra, = -ina^ |a;2. - ^ [B^G%, - B*''G,,] | , (3.27b) 
it becomes easy to see from the self-evident identities 

d^'dT"^^ = (3.28) 
that both real-valued one-forms ^""^i, 

^ + ^ [B^'G*^, - B*^G,,] (3.29a) 

''X = k2u-j^ [B''G\,, - (3.29b) 
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obey a true continuity equation, i.e. 



d^^'^X^, = 



(3.30) 



Therefore they can serve to define the desired normahzation conditions as 

= 1 (3.31) 

is) 

where the hypersurface (S) may be chosen arbitrarily. But clearly, for the present 
stationary field configurations one will choose a time slice {t = const) as the 
hypersurface (5*) which cuts the general form ()3.31|1 down to 

(3.32) 

(3.33a) 
(3.33b) 



rf¥ ^%{r) = 1 
with the modified charge densities ^""^loir) being given by 



B*{r)-X{r) - B{f)-X*{r) 
B*{r) ■ X{r) - B{r) ■ X*{r] 



^X{r) = (2)A;o(f) - 



The second remark refers to the emergence of the gauge potentials 
{*^'*''y4o(r), i?o(r); 74a(r), i?(r)} in the mass eigenvalue equations. In order that this 
system of equations be closed, one has to add the field equations for those poten- 



X{t) ^ {X^{f)} = {Go,im 



221. Thus 



tials, whose general form, however, has been presented already in ref. 
for the present purpose it may be sufficient to merely quote their linearized form 

A («Uo(f) = Anas ^"^Hr) (3.34a) 
A 5o(r) = -A-KUs hl{r) (3.34b) 
Ai;(f) =47ra,4(r1 (3.34c) 



AB{r) 



-Ana^ h*(r) 



(3.34d) 



This simplified form is sufficient for a rough estimate of the atomic energy levels 
(or frequency of spectral lines, resp.). It corresponds to the use of the Coulomb 
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and Breit interactions by the conventional methods in the hterature 
However for our numerical caculations (see below) we will rely upon the exact 
Poisson equations which are non-linear on account of the non-abelian character of 
the two-particle theory. 



D. Energy Functional 

In order to test the quality of the RST predictions in atomic physics, it is not 
sufficient to solve the mass eigenvalue system because the mass eigenvalues MqC^ 
do in general not coincide with the energy Et carried by the field configuration. 
The question of ener gy f unctional has been considered in great detail in some 
preceding papers 13| so that it may suffice here to simply quote the main 
result: 

2 

Et = Y,^-- ^ac^ - ^4'"^ - ^4"*'^ • (3-35) 

a=l 

The general structure of this result is very plausible because it says that the total 
energy Et of any stationary RST configuration is the "corrected" sum of mass- 
energies MqC^. Indeed the necessity for the emergence of the electromagnetic 
(AE',^^™'^) and exchange (Ae!^^^) correction terms arises from the fact that each 
mass energy M^c^ (a = 1, 2) contains already the interaction energy with the 
other particle so that this interaction energy becomes then counted twice when 
one simply forms the sum of mass energies! This is the reason why the mutual 
interaction energy AEt 

AEt = AE^t""^ + AE^T^'^ (3.36) 

must be resubtracted from the sum of mass energies. Since both contributions 
AE^"^^ and Ae!^^^ themselves consist of two terms (i.e. " electric" (e; h) and 
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"magnetic" (m; g)), 

a4""^ = a4'^ + a4"^ (3.37a) 
AE?^) = + AE^?^ , (3.37b) 

the total correction energy AE'y ()3.36j) is built up by four terms: 

AEt = AE^"^ + AeJ"^ + AE^!*^ + ^E^^^ , (3.38) 

which will now be briefly explained. 

The energy correction AE^^ of the electrostatic type is the difference of (the 
sum of) mass equivalents M^^^ and of the electrostatic gauge field energy E^^ ^ i.e. 



AeP =Y.z,- M^^c^ - Eg) , (3.39) 

a=l 



where the electrostatic mass-energy is defined through 



zi ■ m['^c^ = -he j Sr ('Uo(f) ■ ^%a{f) (3.40a) 
Z2 ■ M!f^c^ = -he j £f (^Uo(r) ■ (2)A;o(f) (3.40b) 



Furthermore the electrostatic gauge field energy E^^ is given by 

.(e) 



E'^' = I d'f Ei{f) ■ E2{r) (3.41) 



Similarly the exchange correction energy of the electric type Ae!^'' looks as follows 

2 

a4^) = 5^ ■ Mf ) + 4^*) , (3.42) 

a=l 
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(3.43b) 



(3.43a) 




He 



I 



•r X*{f)-X{f) . 



(3.44) 



47rQ!s 



Obviously both energy corrections of the electric type AEf: and AE\, are of a 
very similar structure but they differ essentially from their magnetic counterparts. 

Indeed, the point with the magnetic interactions is here that the correction 
terms AE^""^ and AE^?^ are built up by the corresponding gauge field energies 
alone, i.e. we have for the magnetostatic case 



One can easily show that this somewhat different structure of the electric and 
magnetic contributions is generated by the combination of the principle of mini- 
mal coupling with the Lorentz covariance of the theory. One can also demonstrate 
that it is just this specific relationship between the electric and magnetic effects 
which brings the RST predictions close to the experimental data. 

In order to see more clearly the specific way in which the cooperation of min- 
imal coupling and Lorentz covariance generate these differences of the magnetic 
and electric energy contributions it may be sufficient to demonstrate this for the 




(3.45) 



and similarly for the exchange case 




(3.46) 
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linear approximation of the RST-Maxwell equations where also the linear Poisson 
equations ()3.34ap - ()3.34d|) do hold. Actually, for this simplified situation, the elec- 
trostatic gauge field energy E^^ can be converted to the corresponding mass-energy 
m'^^c^ (j3.4()a|l - ()3.4nb)l through a simple integration by parts 

2 

2 



4^^^ \ Y.^a-M^:\\ (3.47) 



a=l 



SO that the electrostatic energy correction /S.Ej^ fl3.39p adopts just the value of 
this field energy 

^^T^ \ t^Za-Mi^^c'^E^^K (3.48) 



2 

a=l 



On the other hand, the magnetic gauge field energy can also be con- 

verted to its mass-energy equivalent Mi™''c^, again by simply integrating by parts: 



Am) 



^« ■ 2 

a=l 



\ X^4■Mr)c^ (3.49) 
with 

zi ■ M^^^c^ j ^^^^ ' ^^(^ (3.50a) 

Z2 ■ Mt^c^ j ^2(^) ■ Mr) ■ (3.50b) 

Consequently, the associated magnetostatic energy correction AE^^^ (HTiHI) be- 
comes 

AElp^ ^ -J2za- Mlr^c^ = -E^^ . (3.51) 

a=l 

Evidently the difference in sign between the electrostatic case fl3.48p and the 
present magnetostatic case ()3.5H) is just necessary in order that the total electro- 
magnetic correction AE'^^'"^ (I3.37a|l is built up by the Lorentz invariant product 
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of the four-potentials and currents fca^ 

a4^™) ^ -i^c J d^f ■ (2)^0 -h-A2 + ■ (^Uo - k2 ■ M (3.52) 

= -hic j (fr{^%^{f) ■ (2U^(f) + (%(f) ■ (i)A'^(f)} . 

A similar argument does apply also to the exchange corrections of the electric 
type (AE'^^^) and of the magnetic type (Aii^^^). Namely, one finds again that 
the gauge field energy i?^^ ()3.44|) equals the electric exchange mass in the linear 
approximation 

-\ j^z^.Mi^)c\ (3.53) 

a=l 

SO that the corresponding exchange corrections Ae!^^ ()3.42|) becomes half the value 
hereof 

^4'^ li^^^- ^a^^' = -E^ ■ (3.54) 

a=l 

For the magnetic analogue E^^"^ ()3.46j) one finds in the linear approximation 



4'^ ^ ^ E ^'^ ■ = ^4^ , (3-55) 

a=l 

with the magnetic exchange masses being defined through 

zi ■ Mi^^c^ j ' (3.56a) 

Z2 ■ M^^'^c^ j ^*^^ ■h*ir) . (3.56b) 

With this result, the exchange correction Ae!^^^ ()3.37b|) becomes 

a=l 

= -hic j d^r^B^ir) ■ ho{r) - B{r) ■ h{r) + 5*(f) ■ hl{r) - B*{f) ■ K*(f)} 

= -]^hc j d'f {B^h^ + B*^h*^} , (3.57) 
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which displays an analogous Lorentz invariant structure as its electromagnetic 
counterpart AE^f""^ (ITK^ . 

However as pleasant as the emergence of such a Lorentz invariant structure may 
appear, one should recall that this result is only an approximation. However for 
a comparison of the RST predictions with the experimental data one is interested 
in the highest numerical precision available, and therefore we will rely for our 
calculations below on the original definitions of the energy corrections Ae!^^ dSSUl, 
AeP (HO^ . AeJ™^ (IS33), and AE^^ ^M)- In this way the non-linear and non- 
abelian effects will be fully included in our numerical results (table I and fig.4). 

E. Mass Functional 

Further insight into the nature of the energy functional Ej- ()3.35p is gained by 
a closer inspection of the mass eigenvalues Ma which form the main contribution 
to the energy E^- Especially, one expects that they should contain also the rest 
mass energy (2Mc2) of both particles. 

The desired mass functionals are easily obtained by multiplying both sides 
of the eigenvalue equations for both particles (i.e. equations ()H.2H|l - (j3.25b|l ) by 
the complex conjugate wave amplitudes ^""^R^, ^""^S^ (resp.) and integrating over. 
Thereby the desired mass functionals appear in the following form 

Za-MaC^ = ZlyMc^ + T^^^ + Za-iMi";:^ + Mjf'> + M^"^) + Mf ) + M^^^) . (3.58) 

Here, the first three terms refer to the one-particle contributions and the remain- 
ing mass terms contain the two-particle interaction energies. Observe that the 
latter kind of contributions has been defined already in equations ()3.40a|l - (j3.40b|l . 
(I3.43a,jl -()3.43bjl. (j3.5()a,jl -(;3.5nb;i and (|3.5(ia,jl -()3.56bjl. However, the one-particle 
contributions require now some additional explanation. But recall that each mass 
eigenvalue M^c^ (a = 1,2) ()3.58|) contains the full interaction energy with the 
other particle so that this would be counted twice if one omitted the correction 
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terms Ai?^"'' for the sum of mass eigenvalues ()3.35|) ! 

After the two-particle contributions are clarified now in detail, it is instructive 
to inspect also the one-particle terms. First observe here that the rest-mass energy 
(Mc^) appears here in combination with the renormalization constants Z^^^ being 
defined by 

Zl) = j c?V^,(f)7/;,(f) . (3.59) 

These constants cannot be put to unity (i.e. Z'^^^ =^ 1), otherwise the non- 
relativistic limit would become incorrect (see the discussion of this point in ref . ) . 
Actually, the renormalization constants ^(a)^ are found to depend upon the ki- 
netic energy T(a) of the particles via 

Next, turn to the single-particle interaction energy Eet^ which specifies the elec- 
trostatic interaction energy of the two-particle system as a whole with the external 
source (e.g. nucleus) 

Eif = -he I d'f (^^■Uo(r-) ■ jo(r) = ^-^te ■ (3-61) 

a=l 

Here the (electrostatic) external source is characterized by the time-component 
(^(^^)Ao) of its four-potential ^'^^^Af^, cf. ()2.1|1 . and the total charge density jo{r) is 
the sum of both single-particle contributions ^")/co(r), cf. fl2.17|l . 

m = ^%(r) + (%(r-') , (3.62) 

with the single-particle currents ^^^ko{r) being defined in terms of the wave function 
ipair) through equations ()3.18aj) - ()3.18bj) . 

Surely, the most intricate one-particle contribution to the mass functional M^c^ 
fl3.58|) is the kinetic energy T(a). With respect to the spherical polar coordinates 
(r, 'd, if) , this contribution is the sum of three terms 

r(a) = '■"■^Tr + ('^>r^ + (3.63) 

45 



which refer to the motion in the radial, longitudinal, and azimuthal directions, 
resp. In terms of the ansatz functions ^°'^R±,^"''^S± ()3.9p - ()3.10p . the radial part 
reads 

271 J or or r 

(3.64) 

The longitudinal part ^"-*r^ turns out as the most complicated one and looks as 
follows: 



he fd^r, ^"^Rl d 



271 J r ei-i^ny/^^ d'd 



And finally, the azimuthal part '■"-'Tl^ of the kinetic energy reads in terms of the 
azimuthal current '^"^/c^ ()3.2()c|l : 



2 J r siwd 27i J rsin-i? 1 9^9 9<y5 

dip dip ) 

(3.66) 

With the eigenvalue equations and the explicit form of the energy functional Et 
being at hand now, one becomes able to test the potentiality of RST in the field 
of the two-electron ions. More concretely, one will first solve the mass-eigenvalue 
equations (j3.23j) - (|3.25bj) . etc., together with the Poisson equations ()3.34a|l - (j3.34d|l : 
then one will normalize the solutions according to the prescription ()3.32p . and 
finally one will compute the desired atomic energy level Et by considering the 
value of the energy functional Et fl3.35|) upon the obtained solutions. However, 
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in order to compare these results rigorously to both the experimental data |17[ 
and to the results of other theoretical approaches (e.g. 1/Z-expansion method 
all-order technique in many-body perturbation theory we will treat the 

original non-linear and non-abelian eigenvalue problem whose linearized version 
()H.H4aj) - (j3.34dj) we presented here only for the sake of brevity. 



IV. TRIPLET STRUCTURE 



If RST is equipped with physical meaning, this theory must be able to demon- 
strate the emergence of two kinds of helium: namely (i) ortho-helium (parallel 
spins) with a triplet structure of its energy levels, and (ii) para- helium (anti-parallel 
spins) with a singlet structure (see also fig.3). This means that the totality of so- 
lutions of the RST eigenvalue system ()3.23|l - (j3.25b|l must be shown to subdivide 
into two subsets, the triplet solutions (J = 1) and the singlet solutions (J = 0), 
where however the one triplet member ^Si splits itself up into the three sublevels 
with conventional notation IS* = 1; 5"^ = 0, ±1) being degenerated when no exter- 
nal magnetic field ifgx is present. Thus the first task is to identify that subset 
of solutions to the RST eigenvalue system which corresponds to the conventional 
triplet states \S = 1; Sz = 0,±1). 



A. Triplet fields \S = 1; 5^ = ±1) 

It should be evident that, when the total magnetic quantum number Sz 
(= ■s^(i) + Sz(2)) is extremal (i.e. Sz = ±1), then both single-particle spins Sz(a) 
should point either to the positive or negative z-direction (5^(1) = 5^(2) = ±|). 
Accordingly, restricting ourselves to the positive z-direction, we reparametrize for 
this case the eight ansatz functions ^"'^R±, ^'^^S± in our general ansatz ()3.9|) - ()3.10|) in 
terms of only four real-valued wave amplitudes -R±(t', "(9), 5'±(r, t?) in the following 
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way: 





= -(1)5+ 


2 




(4.1a) 




= -(i)5„ 


- , 
2 


■R- 


(4.1b) 




= -(2)5+ 


- _ ^ , 
2 




(4.1c) 




= -(2)5_ 


- . 
2 


■5_ . 


(4.1d) 



Thus the Pauh-spinors (")0±(r) ()H.8|) - ();il()j) become 



1 1 



(%,(f=^ ^ i?+-C'o' (4-2a) 



1 1 



(1) 0„(f) ^ -C'l' (4.2b) 

(2) 0+ (f) ^ 5+-Cy (4.2c) 



1 1 



(2)0_(f) ^ 5^.(5. 2 _ (4_2d) 

Once the wave functions ipai'f^) have thus been specified, one can turn to the 
associated ortho-currents. From our equation ()3.19|) one easily deduces the charge 
densities ^'^^ko{r) as 

^% = ^{R^ + R^} (4.3a) 
47r 

= ^{S^ + S^} . (4.3b) 
Air 

Similarly, one finds from the component equations ()3.20a|) - ()3.20cp that both cur- 
rents ka{r) ()3.17a|) encircle the z-axis in the same direction, i.e. 

ki = — R+R^ sin-d ■ = ki{r,d) sm.'d ■ (4.4a) 

k2 = — S+S^ sin'd ■ fL = k2{r, {}) sin-i? ■ , (4.4b) 
zn 

see the introductory figure. And finally the exchange density hQ^r) ()3.21|) becomes 

ho = ■^{R+S+ + R.S.} , (4.5) 
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with the corresponding exchange current h{f) ()3.22aj) - ()3.22c|) being found as 

h = —{R+S_ - R.S+} e.r + —{R+S^ + R^S+} 4 (4.6) 
4:71 An 

4= i r]{r, ^) Cr + /i(r, 'd) sin {} . 

Moreover, since the gauge potentials are coupled to the currents, see the linear 
approximation ()3.34a|) - ()3.34d|) hereof, the functional form of the potentials will 
be qualitatively the same as for the currents, i.e. we assume that the radial and 
longitudinal components of the gauge fields vanish i}-^'^ = '^"M^ = = 0) and 
are then left with 

Aa{r) = r sin Aa{r,{}) ■ e<^ (4.7a) 
B{f) =i(3{r,i9) ■ e; + rsini9 5(r,^9) ■ . (4.7b) 

Here we have assumed that all triplet objects depend exclusively upon the radial 
variable r and spherical polar angle so that the triplet configuration is SO (2) 
symmetric around the z-axis. 

B. Mass Eigenvalue Equations (S = l,Sz = +1) 

With the general functional form of all static triplet fields being specified, one 
can introduce these now into the general mass-eig envalue system (pr^ - ()3.2fid;i 
in order to obtain the corresponding eigenvalue system for the triplet states. The 
point here is that, through the triplet ansatz ()4.1a|) - ()4.1d|) . the set of eight complex- 
valued wave amplitudes ^"^R^^r, ^9, </?), (")^±(r, ^9, ^) becomes reduced to only four 
real- valued wave amplitudes -R±(r, '^), S±{r, i!}). Consequently, there arises for each 
particle an ordinary Dirac equation for the radial motion, which is however com- 
plemented by a separate equation for the longitudinal degree of freedom (described 
by the polar angle ??). Thus the radial equation for the positive-energy amplitude 
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of the first particle (a = 1) is found as 



Or 



-rsin^^[A2- R+ + B ■ S+] = - ^^^^ c-R. , (4.^ 



which is to be complemented by its longitudinal counterpart 
OR. 



r^sm^cosd[A2- R+ + B ■ S+\ . (4.9) 

Of course the same procedure does apply also to the negative-energy amplitude 
-R_, which yields the radial equation as 



dr r 



-r sin^ ^[A2-R- + B- S.] = ^ c ■ R+ , (4.10) 



and its longitudinal counterpart as 
OR- 



--r^sm^cos^\A2- R- + B ■ SJ\ . (4.11) 
ov 

The analogous equations for the second particle (a = 2) can easily be found from 
here by simply applying the particle permutation operation flH.2(iaj) - (j3.26dj) and 
are then found to look as follows for the positive-energy amplitude 5"+: 



dr 



+ [(^"Uo + ('Uo] -S^ + Bo-R.+p-R^ 



-r sin^ ^\Ai-S+ + B- R+] = _Ml±^c ■ 5_ (4.12a) 

n 

(99 

—^ = r^sm^cos^[A^-S+ + B-R+] , (4.12b) 
ov 

and similarly for the negative-energy amplitude S- 

- p""Uo + ('Uo] ■ S+- Bq- R+- (3- R^ 



dr r 



+r sin^ ^[Ai-S_ + B- R_] = Ml^Lc . S+ (4. 13a) 

dS 

—^ = r^sm^cos^[Ai-S^ + B-R] . (4.13b) 
ov 
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As was mentioned already in connection with the general form of the mass 
eigenvalue system (see end of Sect. 3c), the preceding eigenvalue equations ()4.8p - 
()4. 13bj) must be complemented by the Poisson equations for the gauge potentials 
^"^y4o(r, -i?), y4a(r, "(9), B{r,'d) and P{r,'d). For the high-precision calculations in 
atomic physics, the linear approximations ()3.34aj) - (j3.34d|l are not sufficient and 
one has to evoke the exact non-linear form of the Poisson equations, see ref . . 



C. Lowest-Order Approximation 

It should not come as a surprise that it is very hard (if not impossible) to find 
analytic solutions of the preceding mass-eigenvalue problem for the ortho-system. 
Therefore one must be satisfied with discussing the problem of degeneration of 
the triplet states IS* = l;Sz = 0, ±1) in an approximative manner. This may 
be done by considering the value of the energy functional Et fl3.35j) upon cer- 
tain approximate solutions of the present eigenvalue problem, both for the present 
states IS* = 1; Sz = ±1) and for the state 15* = 1; = 0) which is to be spec- 
ified hereafter. In this way one can then decide, at least in such a lowest-order 
approximation, whether or not all three states \S = l]Sz = 0, ±1) are actually 
degenerated in RST. 

Our choice of lowest-order approximation simply consists in neglecting all the 
interelectronic interactions, i.e. we put to zero the gauge potentials ^"'^Ao{r,'i!}), 
Aa{r, ■&), BqIt, ■&), (3{r, and B{r). This results in a great simplification of the pre- 
ceding eigenvalue system ()4.8p - ()4.13b|) : namely the angular equations ()4.9|) . ()4.11|) . 
fl4.12bjl and (j4.13b|l merely imply that the two-particle wave amplitudes R±, S± 
become independent of the polar angle (i.e. R±{r, S±{r, R±{r), S±{r), 

and furthermore the radial equations ()4.8|1 . ()4.1()j) . ()4.12a|l and ()4.13aj) are reduced 
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to the following single-particle form 

!^ + (.^U„.fi_ = -^£L±i£e.fl_ (4.14a) 

dr n 

^ + - <"U„ . fl, = fl, (414b) 

dr r h 

^ + (-Uo-^„ = -^^^c.^„ (4.14c) 
dr n 

— — + -S-- ^^""'Aq ■ S+ = c ■ S+ . (4.14d) 

dr r n 

The solutions for these two decoupled one-particle configurations are well-known 

in the literature (see, e.g. ref. |2]J) and may be parametrized by the principal 
quantum numbers (a = 1,2) of both particles. For instance, the one-particle 
mass eigenvalues Ma are given by 

M 

Ma = . . . ..... =^ ■ (4.15) 



{ria = 1,2,3,... 



Thus one can determine the approximate two-particle energy spectrum 
{ET{ni,n2), say) by taking the value of the exact energy functional Et ()3.35|1 upon 
the approximate solutions {R'^,S^} of the reduced eigenvalue system ()4.14ap - 
()4.14d|l . yielding the desired energy spectrum in a lowest-order approximation. 

Naturally one expects that, for such a rough estimate, it will also be sufficient 
to resort to the linear approximation of the Poisson equations for the gauge po- 
tentials, as specified by equations ()3.34a|l - l3.34djl . Using the standard boundary 
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conditions, the formal solutions of these linear Poisson equations are given by 



WAo(r) = -a, ■ / d^f' -^^^ (4.16a) 

la(r) = -as - [ (4.16b) 
B,{f)=a,.Jd\-'^^^ (4.16c) 

5(f) = a.- jd'r'y^^. (4.16d) 

Here one concludes from equation ()4.5|1 that the exchange density /iq is real- 
valued (/lo = /^o) which then also holds for the "electric" exchange potential BQ{r) 
fl4.16c|l . With this lowest-order arrangement, the electrostatic energy correction 
A^;^""^ becomes 




e 



2 




rf¥(i¥ ' ' ' ' , (4.17) 



^ {Rl{r) + Rl{r)) ■ {Sl{r') + ~Sl{r')) 



IQtt .1 .1 \r — f 

and similarly for the electric exchange correction Ae!^^ fj3.54|) 




e' f f^s^^3^> iR+ir)S+ir) + /2„(r)^4r))(/?+(r')^+(r') + i?4r')^4r')) 



167r2 / / \f—f 



(4.18) 

Observe here that the non-relativistic limit (where the "negative-energy" compo- 
nents -R_, S*. are neglected against their positive-energy counterparts -R+, 5*+) of 
this exchange correction AE^^ just coincides with the conventional exchange in- 
tegral AE^"'"^ (jS3I)- The wave amplitude -R+(r) (S'+(r)) corresponds to the first 
(second) particle's wave function ^/'i(f) {ip2{'f^), up to the angular normalization 
factor of (47r)~^/^ ()4.3ap - ()4.3b|) . Thus, as far as this lowest-order approximation is 
concerned, RST is found to be in agreement with the Hartree-Fock approach, see 
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the discussion of the level shift of ortho-helium caused by the exchange interactions 
in fig. 3. 

Both energy corrections A^;^*"^ fPTTTjl and /^E^ (p38|l of the "electric" type 
are much greater than their "magnetic" counterparts AE'^™'' ()H.5H1 and Ae!^^ 
Therefore it makes sense to define an "electric" degeneracy for the triplet 
states 15* = 1; Sz = 0, ±1 >, namely through the requirement that their energy be 
identical only upon the kinetic and electric parts {Et, say) of the exact functional 
Ej^'. 

2 2 

Et = E{^(a) ■ + T(„) + 4 ■ Mi-)cn + IY1 ^^(^a^ + ^i'^)^' • (4-19) 

a=l a=l 

In contrast to this, the complete (but linearized) energy functional would contain 
also the "magnetic" terms, i.e. 

2 

a=l 
1 ^ 

+ ^a{Mi^) + M(-) + M f) + M(^) y , (4.20) 

a=l 

where the renormalization constants Zi^a) ()3.59j) read in terms of the wave ampli- 
tudes R±, S± of the ortho-system (p3a|l - ()4.2bjl 

/C2 _ C2 
dV2±^. (4.21b) 

Finally, the kinetic energies T{^a) ()3.63|) of both particles have to be specified 
in terms of the wave amplitudes S± in order that the total energy E^ ()4.19p 
is completely expressed in terms of wave amplitudes and gauge potentials. Thus 
introducing the triplet ansatz ()4.1a|l - (j4.1djl into the kinetic energies ()3.64|1 - ()3.66|) 
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immediately yields for the first particle (a = 1) 




(4.22b) 



(4.22a) 



(4.22c) 



and similarly for the second particle (a = 2) 




(4.23b) 



(4.23a) 



(4.23c) 



As expected, the "orbital energy" To vanishes for both particles 



(4.24) 



which is a consequence of the vanishing orbital angular momentum of the presently 
considered triplet states ^Si. 

Summarizing, the RST configurations of definite spin direction 15*= 1,5'^ = 
±1 > are well-defined now in their lowest-order approximation, namely through 
their eigenvalue equations ()4.14ap - ()4.15p . potentials ()4.16a|) - ()4.16d|) and energy 
functional Et ()4.19j) . Therefore the question of their degeneracy with the states 
IS" = 1, 5*2 = > of indefinite spin direction can be studied by explicitely working 
out, as the next task, all the RST fields (especially the corresponding energy 
functional) due to those curious states which have spin quantum number S = 1 but 
vanishing z-component 5*^ = 0. Since we are satisfied here with a demonstration of 
the degeneracy only in the lowest-order approximation, we can restrict ourselves 
to showing that the approximate mass eigenvalue equations for the exotic states 
5*2 = do agree with the present ones ()4.14a|l - ()4.15|l . because in this case the 
kinetic {T(^a)) and interaction energies {Ma"^c'^) entering the energy functional Et 
fl4.19p are identical. 
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D. Indefinite Spin Direction 15" = l;Sz = 0) 

The proposal for the exotic configurations with indefinite spin direction consists 
in disposing of the wave amphtudes ^"'•'-R±, ^"'^S±, due to the general ansatz 
flH.l()|l . in the following way: 





= = (^^5+ = = 


(4.25a) 






(4.25b) 






(4.25c) 






(4.25d) 






(4.25e) 



Thus the Pauli spinors ^°■^±{r) ()3.9|) become for the exotic states 

(1) (/,_^(f^ = e*^/2 ^^^^^ ^) . ^(+) (4_26a) 

(f) = -z (^^ ^) . ^(+) (4.26b) 

(2) 0^ (f) = e*^/2 5+(r, i9) ■ (4.26c) 
(2)0_ (f) = _i ^_ (^^ ^) . ^(+) . (4.26d) 

Here one assumes that the four wave amplitudes R±{r,-&), S±{r,'d) are unique, 
complex-valued functions over three-space so that the corresponding Pauli-spinors 
become ambiguous with respect to their sign: 

('^Vilr, + 2n) = -('^Vilr, ^, , (4.27) 

which is frequently the case for spinors, because they transform under SU{2) in 
place of 50(3). 

However the corresponding physical densities (as " observables" of the theory) 
must be real-valued and unique. Actually, inserting the present triplet ansatz 
f)4.25a|) - ()4.25e|) into the electromagnetic currents ()3.19j) - ()3.20c|) yields first for the 
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charge densities of both particles the unique result 

= — (4.28a) 

(''>ko = ^+^+ + ^-^- . (4.28b) 

Similarly, the electromagnetic current becomes for the first particle (a = 1) 

(1)^;, = i -+ ^ 4.29a 

27r ^ ^ 

/IN, it. It- + it+lt , , , 

^^^^ = ± ^ ^ , (4.29b) 

and for the second particle (a = 2) 

= , ^+^- ^+^- (4.30a) 
= - + J' ^ " , (4.30b) 
whereas both currents ka{r) have vanishing azimuthal component (a = 1,2): 

= . (4.31) 

Thus the currents ka{r) ()3.17a|) of the exotic states do not encircle the z-axis (as 
the symmetry axis of the field configuration), but they build up a rotational flow 
in the two-planes containing the z-axis, see the introductory figure. A similar 
geometric pattern does emerge also for the exchange current h{f) (j3.17bj) which is 
found to have vanishing azimuthal component ()3.22c|l 

K = ^ , (4.32) 

whereas the radial and longitudinal components (j3.22aj) - (|3.22bjl are obtained as 

hr = ^ {RIS- - R*_S+] (4.33a) 
2tx 

= -— {RIS- + R*^S+] . (4.33b) 

271 
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Finally, the exchange density ()3.2H) is found to be of the following form 

75* C I D* C 

rt I 0+ + rt_o_ 
ht) = — . 



(4.34) 



Concerning the functional form of the associated gauge fields, one will adapt this 
again to the geometric pattern of the currents. Thus one puts to zero the azimuthal 
components of the magnetic vector potentials {^""^A^ = = 0) which reduces their 
general forms ()3.17cj) - (j3.17dj) to 



B(f) = BrCr + B^e^ , 



(4.35a) 
(4.35b) 



with the components ^""^Ar, ^°'^A^, B^-, B^ being functions exclusively of r and t?. 
Clearly these components must now enter the mass eigenvalue equations for the 
exotic spin states {Sz = 0) which are to be deduced again from the general form 
fl3.23|) - (j3.26d|l by inserting there the present special form fl4.25aj) - ()4.25e|l of the 
wave amplitudes. In this way, the mass eigenvalue equations for the first particle 
(a = 1) appear in the following form: 

Ml + M 



dr re*''/Vsint9<9?9 L 
+BoS. + + z ^^^Ar] ■R+ + [B^ + z Br] ■ S. 

dR^ 2„ i d 



h 



c ■ R- (4.36a) 



^ _l 

dr ^ r ~ re'^/^Vsiu'dd'd 



e'^'/^V^R- - pUo + ('Uo] ■ R, 



-BoS+ - - t (2)AJ .R_-\B^-t BJ ■ S_ = Ml_^c ■ i?^ 

n 



(4.36b) 
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For the sake of completeness, the eigenvalue equations for the second particle may 
also be written down explicitly: 

dS+ i d 



dr re^'^l'^y/sm^O'd 



+BIR_ + + I (^U,] ■ 5+ + [5* + z B;] ■ R+ = _Ml±^c ■ S_ (4.37a) 

/ L 



-S 



dr r re^/'^^/sm.'dd'd 



-B*R+ - [(^U^ - t (^U,] -S.-lB^-t B;] ■ R_ = ^\ ^ c ■ S+ . (4.37b) 

In order to close this system, the Poisson equations for the gauge potentials have 
to be supplied; but since we are satisfied for the moment with the lowest-order 
approximation, one can again resort to the linearized form ()3.34aj) - (j3.34djl . How- 
ever, it is important to mention here that with the present shape of the magnetic 
vector potentials Aa{r) ()4.35a|l and B{f) (j4.35b|) the magnetostatic fields Haif) 

Hi{r) = V X Ai(f) - iB{r) x B*{f) (4.38a) 
#2(r) = V X l2(r) + iB{r) x B*{r) (4.38b) 

are of purely azimuthal character, i.e. 

Ha{^ = ^''^H^ir,^)e^ , (4.39) 

so that these states jS* = 1; 5*^ = > of indefinite spin direction cannot carry a 
magnetic moment, in contrast to the states |S' = 1:5'2 = ±1> (see the discussion 
of this point in ref. [18[)- Therefore only the latter states can provide a handle for 
an external magnetic field in order to shift the energy levels (in opposite directions) 
whereas the energy of the exotic spin state jS* = 1; = > remains unchanged 
(Zeeman effect |l|). 

The fact that both field configurations 15 = 1; 5^ = > and IS* = 0; 5^ = > 
have vanishing ^-component of their total spin S lets expect that their geome- 
tries will share certain common features. In order to see this more clearly, one 
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transforms the triplet wave amplitudes R±, S± in the following way: 



R± = 2 Vsin^ ■ R± (4.40a) 
'S± = e'^V^-S± , (4.40b) 

and additionally one defines complex "magnetic" potentials A^, B-t through 

Aa = ("M^ + i ^"^Ar (4.41a) 
B± = B^±i Br . (4.41b) 

By this arrangement, the present eigenvalue system ()4.Hfiaj) - ()4.37bjl is cast into the 
following form: 



dr r d'd 



Mi + M , 

^ -c-R- (4.42a) 



or r r ov 

Ml- M , 
= —— c ■ R+ 4.42b 

n 

^--i-'S++ P^^^o + ^^Uo] ■ 'S- - Bl . )?„ + Ai ■ ^+ + ■ 
or r ov 

^.Ml±R,.'S_ (4.42c) 

n 

_ Ms -M 

~ h 

Unfortunately it is very difficult to obtain approximate solutions to this system, 
otherwise one could have taken the value of the energy functional upon these 
approximations in order to test the degeneracy hypothesis with respect to the 
whole triplet \S = l;Sz = 0, ±1 >. However, we will consider now the RST 
singlet states IS* = 0;S'^ = >, which obey an eigenvalue system very similar to 
the present one ()4.42a|) - ()4.42d|) for which we are able to present explicit numerical 
solutions. 



c-'S+ . (4.42d) 
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V. SINGLET STRUCTURE 



In order to see clearly both the differences and the common features of ortho- 
and para-helium, it is most instructive to inspect the latter type of helium along 
the lines of arguments being used for the preceding discussion of the triplet states. 
Especially one is interested here in seeing the reason why the energy of the singlet 
states (in lowest order) is not shifted by the exchange interactions in RST (and in 
the Hartree-Fock approach), in contrast to the energy of the triplet states. 



A. Singlet Fields 

It seems self-suggestive that for the RST analogue of the conventional states 
15 = 0; 5*2 = 0) the single-particle spins should be anti-parallel, i.e. one tries for 
the general wave amplitudes ^"^S^ (jTTnil the following ansatz: 

= -WS+ = -^e'^ i?+ (5.1a) 

= -W5„ = -le'^i?_ (5.1b) 

^'K = ^'^S^ = ls^ (5.1c) 

(2)i?_ = (2)5_ = i5_ . (5. Id) 

Thus the Pauli spinors ^"■^±(r) ()H.9|) - (|3.1()j) become 

«0+(r) = i?+-Cy (5.2a) 

(1) 0_(f) = (5.2b) 
(%(r) = 5+-cV' (5.2c) 

(2) 0_(f) = ■ • (5.2d) 

Concerning the corresponding densities, one finds for the charge densities 
^"•^ko^r) ()3.19p formally the same expressions ()4.3ap - ()4.3b|) as for the triplet con- 
figurations IS* = l;^^ = +1). Similarly the first particle's para-current ki{r) is 
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identical to the ortho-current ()4.4a|) whereas the second particle's current k2{'r) 
must necesseraly change sign: 

hir) = S+S- sin-i^eo, . (5.3) 

However the most interesting point refers now to the exchange density hQ{r) ()3.21|) . 
Indeed by inserting there the present para-ansatz ()5.1a|l - (j5.1d|l . one finds that the 
exchange density vanishes identically 

ho{f) = . (5.4) 

Consequently, one puts to zero also the time-component of the exchange potential 

Bo{f) = , (5.5) 

since this component is linked to the exchange density ho^f) via the Poisson equa- 
tion ()3.34bjl . see the solution ()4.1(icj) . But when both objects ho and Bq do vanish, 
the exchange energy Za ■ Ma''' (? ()3.43a|l - (j3.43b|l must vanish, too. However, since 
the latter agrees (in lowest order of approximation) with the exchange correction 
IS.e'^^ ()3.54|) and therefore also with the conventional exchange energy I^Eq fl3.4|) 
(see the discussion below equation ()4.18|) ). one encounters here the origin of the 
missing exchange energy in RST and in the HF approach! 

Though both approaches lead to the same effect of missing exchange energy for 
para-helium, the interpretation is however different: 

(i) The conventional interpretation of quantum mechanics is of statistical nature 
and therefore the missing of exchange energy AEc ()3.4j) is traced back to the non- 
identity of electrons with different z-components (= ±|) of their spins (thus 
escaping the Pauli principle). Such different particles are described by the Hartree 
approach, not the HF approach. As is well-known, the Hartree approximation 
does not take into account the exchange effects [l|, |20| . 

(ii) In RST, the exchange forces are considered as being due to real gauge inter- 
actions (similar to the electromagnetic interactions); and therefore the missing of 
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the ("electric") exchange effect is interpreted as being due to the missing of the 
source for exciting the exchange potential Bq. 

Concerning the functional form of the gauge fields, this will be dictated again 
by the present form of the associated currents. Thus, both electric para-currents 
ki{r) ()4.4aj) and k2{r) will imply again the former shape ()4.7a|l of the vec- 

tor potentials v4a(r), but the ortho-form of the "magnetic" exchange potentials 
B{f) ()4.7b|) cannot be taken over to the present para-case. The reason is that the 
exchange current h{r) must adopt now its para- form as 

h{r) =th{r,^) -W^ir) , (5.6) 

with the complex-valued vector field Wp{r) being defined through 

Wp{f) = -e'^^ie^ + i cos^e^) , (5.7) 

but with the strength of current h{r, being formally the same as for the ortho- 
case (HSJ: 

Mr,^) = ^^-^" + ^--^^ . (5.8) 
47r 

Consequently, with reference to the Poisson equation ()3.34djl . one puts for the 
para-form of the "magnetic" exchange potential B: 

B{r) = irB -Wp , (5.9) 

with a real- valued exchange amplitude B = B{r,d). 

B. Mass Eigenvalue Equations 

After the functional forms of all the para-fields are known, one inserts these into 
the general eigenvalue system ()3.23|) - ()3.26d|) and finds then the mass eigenvalue 
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system for the first particle (a = 1) consisting again of a radial part 
dR+ 



dr 



+ [(^"Uo + ('Uo] ■R.+r[B-S+-A2-R, 



+ rcos^'d[A2- R+ + B ■ S+] = - ^^t^^ c-R 



h 

(5.10a) 



+ - - p^'Uo + (^Uo] -R^+rlA^-R.-B-S. 



dr r 



r cos^ ^[A^-R^+B- SJ^ = ^ ■ R 



h 

(5.10b) 



and of an angular part 



d¥i 

= r'^sindcos'd[A2- R+ + B ■ S+] (5.11a) 

—^ = r'^sui^cos'd[A2-R- + B-SJ\ . (5.11b) 
dv 

Observe here that the angular part ()5.11aj) - ()5.11b|) of the present para-system 
exactly agrees with the angular part ()4.9p and ()4.1H) of the former ortho-system, 
although the exchange potential B^if) looks rather different for both cases, cf. the 
ortho-potential B{f) (j4.7b|l and its para-counterpart ()5.9|l . However, in contrast 
to the angular system, the radial equations do appear rather different, despite the 
fact that the first particle (a = 1) has spin up in both cases; compare the ortho- 
equations ()4.8|) and ()4.10|) to their present para-counterparts (j5.10ap - ()5.10b|) ! 

Finally, the radial part of the second particle's eigenvalue system is found as 
dS. 



dr 



[(^"Uo + (^Uo] -S.+rlA^-S^ + B- R+] 

Mo + M 

+ rcosH[B ■ R+- Ai- S+] = — c-S. 



h 

(5.12a) 



'^-S^ - [(^-'Uo + (^Uo] -S^-rlB-R^+A.-S^ 



dr r 



+ rcos'^[Ai-S.-B-RJ\ = ^\ ^ c-S^ 



h 

(5.12b) 
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and the associated angular part looks as follows 
dS+ 



r'^sm^cos§[B ■ R+- Ai- 



(5.13a) 
(5.13b) 



As was mentioned already in connection with the general eigenvalue equations 
in Sect. 3, such a system must be completed by the Poisson equations for the 
gauge potentials in order to have a closed system. However, for the present singlet 
states ^5*0 we are not satisfied with the linearized form of the Poisson equations 
()3.34a|) - (j3.34djl . but will apply their exact form [2^ for our subsequent numerical 
calculations. 

For both electrostatic potentials ^°'^AQ{r,'d) one finds 

A(r, ^Y^^Ao = 47ra,(i)A;o - 2r\l + cos^ ^)AoB^ (5.14a) 
A(r, 0'^^Ao = 47ra,(% + '2r\l + cos^ i^)AoB'^ (5.14b) 



where A{r,'d) is the 5*0(2) symmetric Laplacean 



A(r,^9) 



Q^2 



2 d 1 



9^2 



cot i} 



d_ 



and the difference of the electrostatic potentials Aq is defined through 

. (Mi-M2)c 



Anff) 



h 



pUo(f)-(2Uo(f)] . 



(5.15) 



(5.16) 



Similarly, the magnetostatic potentials Aa{r, {}) ()4.7a|) must obey the following 
equations 



WON 2 9 cot ^9 9 
A{r,d) + -— + 2 



r dr 



k 2 ^ 

A = 47r«,^ ± {— - — (cos^E^) + 2BYr} 
r sm w ov 



(5.17) 

with the current strengths ka{r,i}) being deduced from their three- vector form 
(jOa|l and (Q as 



ki = —R+R_ 

271 + 



2n 



(5.18a) 
(5.18b) 
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Furthermore, the radial component Yr of the "magnetic" exchange field Y{r) ()3.4fi|l 
reads in terms of the exchange potential B 



Yr = ^7^(cos' ^ ■ 5) + ^'^i/2B (5.19) 

sm V cos V ov ' 



(Ai/2 = Ai- A2) . 

And finally, the Poisson equation for the "magnetic" exchange potential -B(r, §) 
is given by 



cos » a^y r 

-B (A^ + - (2)//^]) . (5.20) 

Here, the exchange current density /i(r, -i?) is given in terms of the wave amplitudes 
by equation ()5.8j) and the objects ^"^iJ^ {ci = 1? 2) are the radial components of the 
magnetostatic fields Ha{r)- 



C. Ortho/Para Level Splitting 

For a discussion of the energy splitting between the present para-configurations 
^Sq and the former ortho-configurations ^Si of Sect. 4, one will resort again to 
the lowest-order approximation. To this end, one neglects all the interparti- 
cle interactions for the para-eigenvalue system ()5.1()a|l - ()5.13b|l and thus finds the 
present system being reduced to the former approximative system fl4.14a|l - (j4.14d|l 
which therefore represents a common approximation to both the ortho- and para- 
configurations. However, both types of solutions are different with respect to their 
energy content Et-, even for the lowest-order approximation Et ()4.19jl . The reason 
is here that the present para-configurations have vanishing exchange density ho{'r), 
cf. ()5.4j) . and therefore the exchange energy AE?^ (jT^ and pTRjl must nec- 
essarily also vanish! Consequently, the para-levels do not receive a shift through 
the exchange interactions, as is the case with the ortho-levels (see fig. 3 and the 
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discussion below equation ()4.18|) ). A further consequence is the fact that the RST 
(and HF) predictions for the energy difference between the singlet states ^5*0 and 
the triplet states ^Si are (in lowest-order approximation) half the value of the pre- 
dictions of the standard theory, which, on the other hand, are in good agreement 
with the experiments; see the discussion of this point in ref. jisj l. 

Does this mean that RST (as the relativistic generalization of the Hartree-Fock 
approach) can be useful only for the computation of the triplet level systems but 
becomes useless when one wants to consider the singlet levels? This question can 
be decided only by concrete calculations of the frequency of spectral lines for the 
transitions within the singlet level system ^5*0 and comparison to the experimental 
data. The point here is that the frequency of a spectral line is computed from an 
energy difference Et^.^.^ — Et^^.^^ due to the total energy Et ()3.35j) of the initial 
{in) and final {fin) states; and it may well be that the influence of the (missing) 
exchange energy AE^\ as part of the total energy Ej-, cancels or at least is weak- 
ened when forming that energy difference. For instance, a closer inspection of the 
" semiclassical" (i.e. u = 0) RST prediction of the ground state (Is^ ^Sq) energy 
Et has revealed that there is a deviation from the predictions of other theoretical 
by not more than 1, 3 eV in the range 2 < z^x < 100, see 
ref.[l6|. Therefore one has to resort now to numerical calculations in order to 
decide the question of precision within the singlet level system. 

VI. NUMERICAL RESULTS 

For a numerical test of RST, one will first consider a sufficiently simple situation 
for which the computational effort remains moderate. In this respect, a consid- 
erable complication of the preceding eigenvalue system is due to the fact that 
the two-particle wave amplitudes R±, S± are functions of both polar coordinates 
{r, i}} which thus yields a two-dimensional eigenvalue problem. Therefore a pos- 



approaches 



10 



n 
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sible simplification can surely be attained by considering isotropic configurations 
where the wave amplitudes {R±,S±} and all the gauge fields {^'^^Aq, Aa, B} de- 
pend exclusively upon the radial variable r. However such isotropic configurations 
are possible as exact solutions only for the singlet states ns^ ^5*0 (Sect.V) where 
the spins of both particles are antiparallel and their principal quantum numbers 
do coincide {rii = 71-2 = n). But in the general case (ni 7^ 71.2), the angular parts 
fl5.11ap - ()5.11b|) and ()5.13a|) - ()5.13b|) of the para-eigenvalue system demonstrate that 
the wave amplitudes must be expected to depend also of the polar angle t?; unless 
one can find solutions for which the right-hand sides of those equations do vanish. 
Obviously this requirement implies the following conditions: 

(6.1a) 

(6.1b) 

(6.1c) 

. (6.1d) 



A2-R+ + B-S+ = 
A2-R^ + B-S^ = 
B ■ R+ - Ai ■ S+ = 
B ■ R_- Ai- S. = 



A. Isotropic Configurations ns^ ^Sq 

Indeed, such isotropic solutions are possible, namely by putting 

(^Uo = = (PUo(r) (6.2a) 

Ai = -A2 = Apir) (6.2b) 

B{r) = Ap{r) (6.2c) 

Wko = (2)A;o = ^P^koir) (6.2d) 

ki = —k^ = h = kp{r) . (6.2e) 

Furthermore both particles are to be described by identical wave amplitudes (-R+ = 
5"+; /2_ = SJ) because one must assume that both eigenvalues do coincide 

Ml = M2 = M" . (6.3) 
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This means that the principal quantum numbers ni,n2 of the conventional no- 
tation must agree {rii = 722) so that one arrives at the RST counterparts of the 
conventional states ns^ "^Sq. These states are the solutions of the following isotropic 
form of the more general system (j5.1()a|l - ()5.10b)l 

^ + r(-Uo + (^'Uol ■ R. + 2rBR^ = - ^"t^ cR. (6.4a) 

+ 2 ^_ _ r(-U„ + (pU J . _ 2rBR^ = i?+ . (6.4b) 

dr r h 

This simple (but exact) eigenvalue system for the isotropic states must be com- 
plemented again by the isotropic Poisson equations for the gauge potentials: 

2d 



^ 2 + -:r ^''^Mr) = 4vra, ^^^oir) (6.5a) 
dr^ r dr J 

(tI + B{r) = 47r a> - 65^(1 - Vi?) , (6.5b) 
\ dr^ r dr J r 6 

which of course is deduced from the original system ()5.1()aj) - ()5.10bj) by means of 
the isotropy conditions ()6.2a|) - ()6.2d|) . 

Naturally, the assumption of isotropy will simplify also the (exact) energy func- 
tional Et ()3.35p . The first point to observe here refers to the normalization con- 
ditions ()3.32j) . Since the "electric" exchange field strength X{r) ()3.33aj) - (j3.33bj) is 
given quite generally in terms of potentials through j22] 

X(r) = -Vfio(r) + i Ao(r) ■ B{f) + i Bo{r) [li(f) - A2(f)] , (6.6) 

this vector field must vanish {X{f) = 0) because both the exchange potential Bq 
and the electrostatic potential difference Aq ()5.16|1 is zero due to the assumption 
of isotropy for the singlet states ns^ "^Sq. Therefore both charge densities ^°'^ko{r) 
and ^^^lo{f) ()3.33a|) - ()3.33b|) do coincide and this puts the normalization parameter 
Za (see below equation ()3.40b|) ) to unity {za = 1) for both particles (a = 1, 2): 

j d^f ^P^ko = j drr^ {Rl + Rl) = l . (6.7) 
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The next point is even more intimately connected with the property of isotropy. 
It refers to the relationship between the energy corrections AE'^^^ and AE^"* of the 
magnetic type. First recall here that on behalf of the vanishing "electric" exchange 
correction Ae!^^ ()H.42j) the exact energy functional Et reduces to 



E' 



T 



E'^ = 2 M"c^ - AeP - AE!f'> - AE't , 



.(9) 



(6.^ 



with the "magnetic" corrections AE!f^^ and Ae'^^ being given by equations ()H.45j) - 
flH.4fi|l . However one can show [3| that for the isotropic states both "magnetic" 



corrections AE'^'^ and AE'l?^ obey the numerical relationship 



AeI^'^ = 2AE^^ , 



-,(m) 



(6.9) 



so that the energy functional for this type of para-states becomes further simplified 
to 

E'^ = 2 M"^ - AeP - 3AE^™) . (6.10) 

Here the magnetostatic correction Ae':^'^ (I3.45|1 reads in terms of the spherically 
symmetric wave amplitudes R±{r) 



AE, 



(m) 



T 



— he 

3 



00 -| 



dr B(r)R<(r)R-(r) 



a 



dr r^ B%r)[l ~r'B{r)] 



s Jo 



(6.11) 

Finally it remains to specify the electrostatic energy correction Ae!^^ for the 
"isotropic" functional -E^ ()6.10|) . Actually, there occurs also a simplification for 



this contribution since the general form of electrostatic fields 



22 



Ei(f) = -V (^Uo(r) 
E2{r) = -V (2Uo(f) + z 



fio(r)B*(r)-S*(f)5(f) 
Bo{f^B*{r)-B*{r)B{r) 



(6.12a) 
(6.12b) 



is reduced for the isotropic configurations {Bq = 0) to the form of simple gradient 
fields 

- -V("Uo(r) . (6.13) 



^a(r) 
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As a consequence, the electrostatic field energy E^^' ()3.4ip equals exactly the elec- 
trostatic mass-energies Mi'^^c^ ()3.40a|) - ()3.40b|) . so that the corresponding energy 
correction Ae!^'' i]'S.'S9\i becomes 

1 ^ r 
a4'^ ^-Y^Za- M^'^c^ = -he dr (^'Uo(r) {Rl{r) + i?!(r)} . (6.14) 

a=l 

This completes the "isotropic" energy functional E^ ()(j.l()|) whose values {E!p\ say) 
upon the solutions i?^^ of the "isotropic" eigenvalue problem ()6.4|) - ()6.5p with mass 
eigenvalues M" constitutes the spectrum for the RST states due to the conventional 
classification ns^ ^Sq. 



B. Relativistic Energy Difference ^Sq / Is^ ^Sq 

Obviously the isotropic form of the RST eigenvalue problem is sufficiently simple 
in order to solve numerically the corresponding eigenvalue equations ()6.4a|) - ()6.5b|) 
and to take the value of the exaci energy functional -E^ ()6.10jl upon these solutions. 
This admits us to calculate, e.g., RST energy differences within the spectrum of 
isotropic solutions and to compare this to both the experimental data and to the 
analogous predictions of other theoretical approaches in the literature. As a simple 
demonstration, we now inspect the energy difference of the states "^Sq (i.e. 
rii = 77,2 = 2) and Is^ -"-So (i.e. ni = n2 = 1), where the latter state is the ground- 
state of the helium-like ions and has been studied already extensively The 
energies to be considered refer to an idealized situation were the nucleus is thought 
to be infinitely heavy and the self-energy effects (and other QED corrections) are 
first neglected (for the inclusion of self-energy in RST, see ref.s jisl Q|)- Thus 

o 

denoting the semiclassical energy difference by A Ei\2, we have 

o o n 

AEi\2 = E 

^ 2s2 iSo 

where the energy functional to be used in RST is that one of the isotropic states 
fj6.10|) . This RST prediction may then be compared to the observational data 
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(6.15) 



{AexpEi\2) as compiled by NIST 



irj l. Such a comparison is instructive because 
these experimental values Aea;p-E'i\2 include the Lamb shift which in the present 
("semi classical") RST approach is first neglected and afterwards taken into account 
by switching on the RST self-interactions. But despite this neglection it is interest- 
ing to inspect also the corresponding semz-c/asszca/ predictions of other theoretical 



approaches, especially the relativistic (l/Z)-expansion method 



Ai/z-E'i\2) 



and the all-order technique in relativistic many-body perturbation theory (MBPT) 

Q (^ AallEl\2). 

The latter approaches are not directly concerned with the special energy dif- 
ference A Ei\2 but rather with the (semiclassical) one-particle ionization 
energies, Ji^^{2) and J^n (2), namely for the situation when the principle quan- 
tum number of one particle is n, and the other particle is in the ground state 
[n = 1). However it is an easy matter to compute the desired energy difference 

o a{n) o(n) 

A Ei\2 from those helium-like {Jtheo (2)) and hydrogen-like {Jtheo (1)) ionization 
energies. Indeed, when one of the two electrons is thrown out from a bound state 

o(n) 

of energy Efheo (2) to infinity, where its energy is just Mc , the remaining bound 
electron has energy Etheo{^)'i ^^us the ionization energy Jt/jgo (2) emerges as 



the difference of the initial and final energies: 



o(n) 
Jtheoi'^) 



n o(n) 

Mc'+ E,,,M) 



o(n) 

EtHeoi^). (6.16) 



o(") 

Here the semi-classical one-particle spectrum Etheo (1) is due to the simplified 
eigenvalue problem ()4.14a|) - ()4.14b|) which yields 

o(n) ^ ^ ~ „ MC^ 

E,,,,(l) ^ M^c' = , == . (6.17) 

1 + 



^n-l+y/l-(zexas)^ 

Now one can write down the equation ()6.16|) for both principle quantum numbers 

o 

n = 1 and n = 2 and then finds the desired energy difference Atheo Ei\2 by 
subtracting both equations from one another: 

o(l) o(2) o(2) o(l) 

Ai,.eo^l\2=Ji,,eo(2)- Jtheoi'^)+ ^i/.eo(l)- Etheoi^) ■ (6-18) 
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Here one can now insert the one-particle energies from the semiclassical spectrum 

o(n) 

()6.17|) and the two-particle ionization energies J^^^g^ (2) for n = 1 and n = 2 



from ref.s jlfll . in order to yield {theo l/Z^all) the energy difference 

o o 

\/zEi\2 of the 1/Z-expansion method and \11E1\2 of the all-order technique in 

o(n) o(n) 

MBPT. Unfortunately the ionization energies J^/^ (2) and J^n (2) are tabulated 
in ref.s jiol . exclusively for the states Is^ "^Sq and ls2s -^Sq but not for the 
presently considerd 2s^ "^Sq states, so that we have to restrict us preliminarily to a 
comparison of our RST results Ajist Ei\2 to the experimental data Aexj^i\2 12\, 
see table I and fig. 4 . 



Table I 

C. Self-Interactions 

One could suppose that the origin of the deviations between the present semi- 
classical RST results Ajist Ei\2 and the experimental values Af,xpEi\2 (table I 
and fig.4) are due to the neglection of the electronic RST self-interactions. This 
however is not true. Indeed, we will readily demonstrate that the inclusion of 
the RST self-interactions is not sufficient in order to shift the semiclassical RST 
predictions Arst Ei\2 (jS.lSj) towards the experimental values AexpEi\2. This is in 
contrast to the situation with the ground state Is^ ^Sq, whose ionization energy 
deviates not more than 2 eV from the experimental values in the whole range 
2 < Zex 83 when the RST self-interactions are adequately included. The RST 
self-energy problem has been treated extensively in ref.s Q] so that it may 
be sufficient to reproduce here only the main results which are relevant for the 
computation of the desired energy difference AjistEi\2: 

ArstEi\2 = Ex — Et . (6.19) 

^ 2^2 iSo ls2iSo 
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Since this looks quite similar to the semiclassical result ()6.15|) without the self- 
interactions, it is important to clearly state now the differences. 

First let us mention that our RST self-interactions enter the theory via the 
fibre metric Kap ()2.2()j) . The general form of this fibre metric is parametrized by 
the "self-interaction parameter" u : Kap = Kap{u)\ and this parameter enters 
then also both the energy eigenvalue system and the energy functional (i.e. 
Et =^ Et{u)). However, for the presently considered isotropic solutions ns^ ^Sq, 
the change caused by a non-trivial value of u is minimal: First, the mass eigenvalue 
system ()6.4ajl - (j6.4bjl remains invariant under a change of u. Next, it is easy to show 
that also the electrostatic Poisson equation ()fi.5aj) remains the same for a non-zero 
value of M. Thus it is only the magneto static Poisson equation (j6.5b|l which changes 
into 

d? Ad 



dr"^ r dr 



+ -— 1 B{r)+QB'^{r) 



2 

-r'^B(r) 



47ra,e-2"-^ . (6.20) 



And finally, the energy functional Ei[, ()6.10|) remains the same but with the mag- 
netic correction energy AE^^ being modified according to 

AEp = — hc\ dr r^B{r)R+{r)R_{r) / dr r^B^r) [l - r^B{r)] \ . 

{Jo Ois Jo J 

(6.21) 

Thus it becomes evident, that the inclusion of the self- interactions implies for the 

2 

isotropic states nothing else than replacing the coupling constant as (= j-) by 
as ■ e~^"; but only for the magnetic interactions whereas the electric interactions 
are unchanged! Since the non-singular solutions B{r) of the magnetic Poisson 
equation ()6.20|) must vanish {B{r) = 0) when the right-hand side becomes zero for 
M — s> oo, cf. ()4.16d|) . the magnetic effects dissapear completely and we are left with 
the purely electrostatic interactions. 

This effect can be easily understood by recalling the specific way, in which 
the self-interactions do appear in RST. First, reconsider the relationship ()2.2()aj) - 
fl2.20b|) between the Maxwell currents j"^ and the RST currents j^^ ■ Obviously, 
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the link between both currents is given by the fibre metric K^p which for the 
present two-particle configurations adopts the following form 15[ for non-trivial 
self- interaction (m 7^ 0): 

/ e^sinhw — e^coshw 
-e" cosh u e" sinh u 











\ 





-2n 



V 



-2n 



/ 



(6.22) 



Therefore the inclusion of the self-interactions (i.e. u ^ 0) generalizes the 
former relationships ()3.18a|) - ()3.18dj) to the following form 



■1 



3 



3 M 

J u 



-e "{coshw ■ ki^ + sinhw ■ k2^} 
-e~"{sinhu ■ fci^ + coshn ■ k2^} 



3 - 3-2^/1* 



— e 



-2ui 



(6.23a) 
(6.23b) 
(6.23c) 
(6.23d) 



Thus, putting here n ^ 00 lets vanish the Maxwellian exchange currents j^^ and 
j^^ so that the exchange potential i?^ must be put to zero in order to satisfy the 



last two Maxwell equations ()2.11d|) - ()2.11e|) . Moreover, the first two transformation 
relations ()6.23a|) - ()6.23b|) degenerate for m — 00 to the following form 



2 



'2(^1^ + ^2m) 



J 



(6.24a) 
(6.24b) 



.e. both Maxwellian charge densities become identical (cf. (lOdj) ): 



3^ 0= f 



(6.25) 



but the Maxwellian three-currents ja = {—j"' k]k = 1, 2, 3} must vanish 



Ji =32 = --{ki{f) + k2{f)) 



(6.26) 
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because for the ns^ ^5*0 states the RST currents ka{r) are antiparallel, cf. ()6.2ep 

ki{r) = -Hr) = kp{r) . (6.27) 

Clearly, when all the Maxwellian three-currents do vanish according to ()6.26|) . 
then all the magnetic potentials Aa{r) ()6.2b|) and B{r) ()6.2c|) are also zero as 
required by the Maxwell equations (j2.7bjl - ()2.7e|l . 

Thus, the result is that for u oo one ends up with a configuration which 
has been called "electrostatic approximation" j^], because all the magnetostatic 
fields have been neglected. Now it turns out that this is an exact solution, 
namely for inclusion of the self-interactions, but with n ^ oo. However, the 
more important point with this type of solution refers here to the magnitude of 
the self-interactions. It is true, for vanishing self-interactions {u = 0) the RST 
predictions for the energy differences AfiSTEi\2 come close to the experimental 
values by less than 0, 3% (fig.4); but when u is increased from zero to infinity, the 
amount of lowering of the energy Et is not sufficient in order to shift the RST 
prediction AjjstEi\2 sufficiently close to the experimental value Aea;p-E'i\2 for small 
and intermediate nuclear charge numbers Zex {zex ^ 50, say). This is the same 
effect as was previously encountered in connection with the ionization energy of 
the helium-like ions [iq] where the value of the self-interaction parameter u had 
to be fixed for the ground-state Is^ ^Sq as u = 0.03052 in order that the RST 
predictions meet with the observational data for the ionization energies in the 
range z^x ^ 50. However, for the smaller values of nuclear charge Zex 50), such 
an acceptable fit of the RST parameter u in order to adapt the RST predictions 
to the experimental data is not possible, neither for the ionization energies of ref. 
nor for the present energy differences AEi\2- 



fig.4 
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In order to present a numerical demonstration of this effect, consider the nu- 
clear charge number = 42, molybdenum (see table I). The two-particle ground- 

o(l) 

state energy for m = is found by numerical integration as Et{2)= 9'''3 577.2 eV 

o(2) 

and the excited state 2s So has energy Et(2)= 1 009 805.2 eV. Thus, for van- 

o(2) 

ishing self-interaction [u = 0), one finds the RST prediction Arst Ei\2=Et(2) 

o(l) 

— Et{2)= 36 228eV. Now switch on the self- interactions and let u adopt its 

o(l) 

ground-state value Ui = 0.03052 which lowers the ground-state energy Et(2) 



-^T(2) 



= 973 575.6 eV, i.e. only 1.6 eV below its value for n = 0! Further- 

Ul 

more, taking the most favourable value ^2 = 00 for the excited state 2s^ ^Sq yields 

= 1 009 804.6 eV, which means a lowering of 



the corresponding energy -ELL 

0(2) 

Et{2) by 0.6 eV. Consequently the desired energy difference remains practically 



(2) 

the same: AfjSTEi\2 =f E) 



-^T(2) 



36 229 eV, when one switches on 



1\2 -r ^j<(2) 

the self-interactions; and thus the semiclassical RST predictions cannot be suffi- 
ciently improved for the excited state 2s^ ^Sq by means of taking into account 
the self-interactions. However, for large charge numbers Zex, the energy shift 

o 

AfisT Ei\2^ AjiSTEi\2 becomes larger (fig. 5) and therefore it becomes possible 
for large enough z^x to select a parameter value U2 such that the RST predic- 
tion AjiSTEi\2 fl6.19|) coincides with the (presently unknown) experimental value 
Aea;p£'i\2- A similar effect occurs with the ground-state ionization energy see 
figs. 4 and 5. 

fig.5 

Of course, a rigorous proof of the claimed coincidence of the RST predictions 
AfjSTEi\2 with the (missing) experimental data Aexi^i\2 of ref. for large Zex 
is not possible, but a rough estimate may give nevertheless a first hint upon the 
correctness of that extrapolating claim. Observe here that the experimental one- 
electron ionization energy j}^^2) '^^ helium ground-state for molybdenum {ze 
42) is jJJ^j2) = 23 791.7 eV, see ref. 



i3 



; and this is roughly the same order of 



77 



magnitude as the present RST energy difference Arst Ei\2 = 36 175.6 eV, see 
Table I. However, not only the orders of magnitude are the same, but also the 
relative deviations A^ist of the RST predictions from the experimental values: 
For the ionization energies, one defines 

7(1) 

Arst- -(I) [b.Zii) 

f;xp(2) 

and then finds for molybdenum {z^x = 42) 



o(l) 

'-^RST 



23 850.0 -23 791.7 

> 0.16% , (6.29) 



=42 23 791.7 

see ref.^^. On the other hand, one may define the analogous deviation for the 
present energy differences as 

o 

°(A2)_ Arst Ei\2 -^exjJ^l\2 , s 

Arst— X (^D.dUj 

^ea:pC/l\2 

which yields for Zex = 42 by means of table I 



o(l\2) 
'-^RST 



36 228,6 - 36 175^6 (6.31) 



Zex=i2 36 175.6 



and this is just the same deviation as for the ionization energies ()6.29|) ! (For a 

o(l) o(l\2) 

comparison of both deviations Arst ^^'^ ^rst functions of charge number Zex, 
see fig.4). Consequently, it seems reasonable to assume that for large Zg^ (^ex ^ 40, 
say) the RST predictions ArstEi\2 can be optimally fitted to the experimental 
valuer A„,£.,, by choosing fa the self-inte.action pa.a,nete^, again =o„.e value 
in the vicinity of its ground-state value Ui = 0.03052, see ref. [16||. Table II presents 
an overview of the corresponding RST predictions ArstEi\2 for U2 = ui = 0.03052 
for the higher nuclear charges {zex ^ 40). Clearly, it must be left to the outcome 
of the future experiments for Zex ^ 40, whether this hypothesis Ui ~ U2 is actually 
realized in nature. 

Table II 
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o{l) 

Et{2) 


o(2) 
Et(2) 


AflST Ei\2 






^rstEi\2 


^expE 




u = 


u = 


= 


11 ~ OO 


11 = OO 


11 — OO 


\17 251 


2 


1021920.0 


1021978.2 


58 


3 


1021920.0 


1029782.8 


58.32 


57.87 


3 


1021800.9 


1021948.4 


147 


5 


1021800.9 


1021948.4 


147.46 


-.- 


4 


1021627.4 


1021904.8 


277 


48 


1021627.4 


1021904.8 


277.44 


-.- 


5 


1021399.4 


1021847.6 


448 


24 


1021399.4 


1021847.6 


448.28 


-.- 


6 


1021116.9 


1021779.9 


659 


97 


1021116.8 


1021776.8 


660.00 


-.- 


7 


1020779.8 


1021692.4 


912 


59 


1020779.7 


1021692.4 


912.65 


-.- 


8 


1023881.3 


1021594.3 


1206 


15 


1020388.0 


1021594.3 


1206.25 


-.- 


9 


1019941.8 


1021482.5 


1540 


67 


1019941.6 


1021482.4 


1540.84 


-.- 


10 


1019440.7 


1021356.9 


1916 


19 


1019440.4 


1021359.9 


1916.54 


-.- 


11 


1018884.8 


1021217.7 


2332 


86 


1018884.4 


1021217.6 


2333.23 


2325.620 


12 


1018274.0 


1021064.7 


2790 


66 


1018273.4 


1021064.6 


2791.16 


2782.300 


13 


1017608.2 


1020897.8 


3289 


63 


1017607.4 


1020897.8 


3290.33 


3280.460 


14 


1016887.2 


1020717.2 


3829 


94 


1016886.4 


1020717.1 


3830.75 


3820.313 


15 


1016111.1 


1020522.6 


4411 


55 


1016120.0 


1020522.6 


4412.61 


4401.250 


16 


1015279.6 


1020314.2 


5034 


62 


1015278.2 


1020314.1 


5035.92 


5022.820 


17 


1014392.6 


1020091.8 


5699 


26 


1014391.0 


1020091.7 


5700.83 


5686.783 


18 


1013449.9 


1019855.4 


6405 


56 


1013448.0 


1019855.4 


6407.40 


6392.000 


20 


1011397.1 


1019340.5 


7943 


37 


1011394.4 


1019340.4 


7945.98 


-.- 


22 


1009119.8 


1018769.0 


9649 


11 


1009116.1 


1018768.8 


9652.64 


9631.588 


23 


1007896.5 


1018461.7 


10565 


25 


1007892.3 


1018161.6 


10569.30 


10546.591 


24 


1006616.5 


1018140.2 


11523 


77 


1006611.7 


1018140.1 


11528.37 


11503.997 


25 


1005279.5 


1017804.3 


12524 


80 


1005274.1 


1017804.2 


12530.03 


12503.930 


26 


1003885.4 


1017453.9 


13568 


52 


1003879.3 


1017453.7 


13574.40 


13546.265 


27 


1002433.9 


1017089.0 


14655 


04 


1002427.1 


1017088.8 


14661.69 


14631.622 


28 


1000924.8 


1016709.4 


15784 


59 


1000917.2 


1016709.2 


15791.99 


15759.754 


29 


999357.8 


1016315.1 


16957 


29 


999349.3 


1016314.9 


16965.59 


16930.909 


30 


997732.6 


1015905.9 


18173 


33 


997723.2 


1015905.7 


18182.56 


18145.086 


36 


986743.2 


1013135.2 


26392 


01 


986726.7 


1013134.8 


26408.10 


26354.080 


40 


978212.8 


1010979.0 


32766 


14 


978189.9 


1010978.4 


32788.48 


-.- 


42 


973577.2 


1009805.2 


36228 


6 


973575.6 


1009804.6 


36229.0 


36175.611 


50 


952478.7 


1004445.1 


51966 


49 


952433.0 


1004444.2 


52011.19 


-.- 


60 


920009.2 


996138.4 


76129 


19 


919928.4 


996136.6 


76208.23 




70 


880058.5 


985819.0 


105760 


47 


879926.4 


985816.0 


105889.63 




80 


831546.6 


973137.6 


141591 


00 


831342.5 


973132.8 


141790.31 




90 


772857.3 


957566.4 


184709 


14 


772554.0 


957558.6 


185004.60 




100 


701428.2 


938258.9 


236830 


70 


700988.6 
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Table I: Comparison of Energy Differences AEi\2 



The energy difference AEi\^2 without self-interaction ()6.15|) f4th column) and with 
self-interaction (|6.19|) (7th column) is compared to the experimental values AexpEi\2 

o(l\2) 

(last column). The relative deviations AjiST are illustrated in fig. 4. In contrast 

to the situation for large Zex, the inclusion of the RST self- interactions cannot improve 
the predictions, because for low and intermediate charge numbers Zex the lowering of 
the energies £'^^(2) RST self-interactions is too small, see fig.5. 
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Fig. 4: Ionization Energies J rst{2) Energy Differences 

L^RST El\2 



o(l) 



CI ■• (6-30) 
A;st:(6.28) 



The relative deviations /S.^sT for the semiclassical ionization energy 

o(l) o(l\2) ___ 

J RST(2) the deviation Aj^sT for the semiclassical energy difference 

Arst Ei\2 are of the same magnitude (< 0, 2%) for large Zex(^ 40). This suggests to 



adopt for the excited state 2s Sq the same value for the self- interaction parameter u 



(— > 0,03052) as for the ground-state Is Sq, see ref. 
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Fig.5: RST Energy e!^^^~^ for the excited state 2s^ 

(2) 

The total energy £'-^^2) ■^ex = 46 (Pd) as a function of the self-interaction 
parameter u decreases too weakly for small and intermediate values of Zex in order 
.... .he e.,e...^fee„ce A«..E,. £13 can be ««ed .o .he expeH.enU. value 
AexpEi\2 of ref.[13] by means of a suitable choice of u. 
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Table II: RST predictions for the energy difference AEi\2 

Adopting the same value u = 0,03052 for both the ground state Is^ and the excited 
state 2s^ -^Sq yields the RST predictions for the energy difference AiistEi\2 (|6.19() of these 
states in the limit of large Zex(^ 40) where observational data are not yet available. 
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